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ABSTRACT 

We derive detailed theoretical models for 1074 nearby stars from the SPOCS (Spectroscopic Prop- 
erties of Cool Stars) Catalog. The California and Carnegie Planet Search has obtained high-quality 
(R ~ 70000 - 90000, S/N ~ 3 - 500) echelle spectra of over 1000 nearby stars taken with the Hamil- 
ton spectrograph at Lick Observatory, the HIRES spectrograph at Keck, and UCLES at the Anglo 
Australian Observatory. A uniform analysis of the high-resolution spectra has yielded precise stellar 
parameters (T e ff, logg, usini, [M/H] and individual elemental abundances for Fe, Ni, Si, Na, and 
Ti), enabling systematic error analyses and accurate theoretical stellar modeling. We have created a 
large database of theoretical stellar evolution tracks using the Yale Stellar Evolution Code (YREC) 
to match the observed parameters of the SPOCS stars. Our very dense grids of evolutionary tracks 
eliminate the need for interpolation between stellar evolutionary tracks and allow precise determina- 
tions of physical stellar parameters (mass, age, radius, size and mass of the convective zone, surface 
gravity, etc.). Combining our stellar models with the observed stellar atmospheric parameters and 
uncertainties, we compute the likelihood for each set of stellar model parameters separated by uniform 
time steps along the stellar evolutionary tracks. The computed likelihoods are used for a Bayesian 
analysis to derive posterior probability distribution functions for the physical stellar parameters of 
interest. We provide a catalog of physical parameters for 1074 stars that are based on a uniform set 
of high quality spectral observations, a uniform spectral reduction procedure, and a uniform set of 
stellar evolutionary models. We explore this catalog for various possible correlations between stellar 
and planetary properties, which may help constrain the formation and dynamical histories of other 
planetary systems. 

Subject headings: planetary systems — stars: fundamental parameters — stars: interiors 



1. INTRODUCTION 

Precise analysis of high-resolution spectra of stellar atmosphere and theoretical calculations of the physical properties 
of low-mass stars are essential for a variety of astronomical problems. Many previous works have focused on stellar 
abundance analyses and stellar age determinations to unders tand the chemical and dynamical history of the Galactic 
disk. Since the benchmark work bv lEdvardsson et al.l (|1993af) . there have b een a number of large observational surveys 
to obtain a true age -meta llicity relationship in the solar neighborhood (|Feltzing et al J 120011 : llbukivama fe Arimotd 
120021 : iNordstrom et al.ll200l . 

There has also been a recently growing interest in the properties of stars with planetary companions. Extensive 
Dopplcr radial-velocity surveys using high-resolution spectroscopy and large transit programs have made tremendous 
progress in the past decade and garnered nearly 200 extrasolar planets to date. The distributio ns of observed planetary 
properties are of gr eat importance for testing theories of planet formation and dynamics ( Eggcnbcr ger et all l2004t 
iMarcv et al.l l2005ah . The extensive spectroscopic observations in searches for planets have also revealed statistical 
differences between the planet-host stars and the normal dwarf stars with no companions. It is now well established 
that the frequency of the giant planets with orbital period less than 3 years is a strong function of the stellar metallicity 
for sol ar-type stars. Consequently, the planet-host stars exhibit a diffe rent metallicity distribution from that of single 
stars (S antos et alj 12001 : iFischer fe Valenti 120051 : ISantos et al.ll2005t ). The stellar atmosphere pollution by planet 
accretion has been propo sed as one of th e possible metal enrichmen t scenarios, and it has been tested by observation s 
and theoretical models (|Gonzalezl Il997t lLaughlin fe Adams! Ii997fc ISandquist et all 119981: iPinsonneault et al.ll2001h . 
The efficiency of metallicity enhancement by planet accretion is still under debate, but the models require an accurate 
knowledge of the mass of the stellar convective zone. Deriving accurate distributions of stellar metallicity, age and 
convective zone size from a large sample of observations is useful for better understanding the origin of the metal-rich 
atmospheres of planet-host stars. 

Theoretical models of planet-host st ars are also r elevan t to the dynamical studies of planetary systems. A motiva- 
tional work for this paper was done bv lFord et al.l (|1999j ) (hereafter FRS99), who provided models for five stars with 
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known "hot Jupiters" to estimate the orbital decay timescales for these systems. Specifically, they constrained the 
models of tidal dissipation and spin-orbit coupling in those systems using the estimated convective envelope masses 
and stellar ages. Age estimates of planetary systems also help constraining long-term perturbations on planets. For 
example, if a planet resides in a wide stellar bi nary system, the planetary orbit may undergo secular evolution, on 
a timescale as large as ~ 1 Gyr or even longer : I lohnan et al.l fl997t iMichtchenko fc Malhotral l2004t iTakeda fc Rasiol 
120051 : iMudrvk fc W u 2006) . Estimates of the age and the physical properties of the host star can thus help constrain 
the dynamical history and the formation channel of the system. 

Despite the recent improvements in high- resolution spectroscopy and data analysis techniques, providing accurate 
stellar parameters is not a straightforward task, even for bright stars in the solar neighborhood. Currently, there is no 
canonical method within data analysis or theoretical modeling to estimate the stellar properties. Various approaches 
have been applied, particularly for stellar age determinations. The observed stellar rotation may be a good age 
indicator, since stars normally slow down as they age. Chromospheric activities measured from the C A II H and K 
absorption lines are the favored rotation measure and therefore often used as a stellar age in dicator (|Wilsonl 119701 : 
iBaliunas et"al]|1995l [l997t iHenrv et all 119971 : iRocha-Pinto fc Macie]|fl998t iHenry et al.l l2000al). The stellar age can 
also be constrained by the surface lithium depletion ( Soderblomlll985 : lBoesgaardlll99l( rthough it needs to be treated 
with extra caution for planet-h ost stars as close-in p lanets may tidally affect the stellar convective envelope and thus 
cause further lit hium depletion (i sraclian et al. 2004j). The merits and challenges of different techniques have been well 
summarized by ISaffe et al] (|2005f ) . Thev have carefully derived the ages of 50 extrasolar planet host stars observed 
from the southern hemisphere and tested several different age determination schemes for comparisons. 

In this paper, we derive various stellar properties by matching the spectroscopically determined surface parameters 
to theoretical stellar evolution models. This approach is similar to the traditional isochrone method. In the isochrone 
method, the observed My and B — V are placed in the HR diagram, then the stellar ages are derived by interpolating the 
observed position between the theoretically computed isochrones (Frwarodfl980l: IVandenberd 1985: Ed vardsson et al 



observ ed position be tween the theoretically computed isochron es (iwarog iy»U; Vandcnbcrs iy»o; hjdvardsso n et al 
1993at iBertelli et al.|[l99l iNg fc Bertellil [1991 lLachaume et alj H9991: llbukivama fc Arimotol l2002t INordstrom etal 



2004T |Pont fc Everll2OOl " lJ0rgensen fc Lindegrenl 120051 ; iKaratas et al.H2005h . Equivalents, stellar ages can be derived 



by interpolating the observations between theoretically computed stellar evolutionary tracks. Isochrone or evolutionary 
track analysis is becoming increasingly accurate relative to other methods such as age - activity relations or age - 
abundance relations, given the availability of advanced high-resolution echelle spectrographs, sophisticated stellar 
evolutionary codes, and increased computational power. Theoretical evolutionary models also have an advantage in 
that one can create a full model of a star, providing not just the stellar age but all the physical parameters, including 
those characterizing the stellar interior. 

Note that the traditional spectral data analysis using theoretical stellar evolutionary tracks (or equivalently, 
isochrones) involves many sources of systematic bias. Here we summarize the main steps in the theoretical mod- 
eling together with the necessary precautions: 

(i) The observed sample — A substantial number of sample s tars ar e required to obtain meaningful stellar param- 
eter distributions. The benchmark work by lEdvardsson et al.l (|1993a| ) provided chemic al abundances of 189 nearby 
field dwarfs. Today, ty pical solar neighborhood surveys include ~1000 to -10000 stars (jlbukivama fc Arimotol[2002t 
INordstrom et a l. 2004) . As the sample size becomes less of a problem, however, the sample selection remains crucial for 
understanding the global statistics of the stellar properties. Unfortunately, any survey is limited by its own selection 
criteria, and in practice any sample has some selection effects. The resultant obs erved statistics need to be analyzed 
very carefully to separate the selection effects from the true stellar properties. INordstrom et al.l (|2004l ) have done 
careful studies on the completeness of their sample stars, in terms of binarity, magnitude, sampling volume and other 
stellar parameters. 

For this work, we have us ed the spectroscopic data from the SPOCS (Spectroscopic Properties of Cool Stars) catalog 
bv lValenti fc Fischerl (|2005l . hereafter VF05). SPOCS consists of high-resolution echelle spectra of over 1000 nearby F-, 
G- and K-type stars obtained through the Keck, Lick and Anglo-Australian Telescope (AAT) planet search programs, 
including the 99 stars with known planetary companions. The sampling criteria for the SPOCS catalog are such that 
the achievable Doppler velocity precision is optimized for planet detections. Thus, the catalog favors stars that are 
bright, chromospherically inactive (or slow-rotating) on the main-sequence or subgiant branch. A subset of the catalog 
can be used for completeness studies in terms of sampling volume or the presence of planetary companions. 

(ii) Spectral data analysis — Currently there is no standardized technique for analyzing high-resolution spectra. 
In the traditional abundance analysis approach, different radiative transfer and stellar atmosphere models produce 
different types of errors. Apart from the choice of model, micro-turbulence needs to be carefully adjusted to match 
the equivalent width of the spectrum, and also macro-turbulence, stellar rotation, and instrumental profile need to 
be taken into account for the spectral line broadening effects. All these factors may induce systematic biases in the 
derived model parameters. 

In contrast to the traditional abundance analysis method, VF05 directly fit the observed spectrum to the synthetic 
spectrum generated by a software package SME (Spectroscopy Made Easy. IValenti fc Piskunovl[l996l ). allowing for a 
self-consistent error analysis. Using SME, VF05 have derived a set of observational stellar parameters my, T c g, [Fe/H] 
and logg for each star, along with precise error estimates. The spectral analysis technique s adopted for the SPOC S 
catalog are summarized in § ??. Fo r more details on the observations and SME pipelines, see IValenti fc Fischerl ((2005) . 

(iii) Model parameter estimate — iPont fc Everl (|2004h have thoroughly analyzed the systematic biases induced by the 
traditional maximum likelihood approach using theoretical isochrones, which simply interpolates the nearest isochrones 
to the observational data point in the HR diagram. The simple isochrone interpolation approach accounts for neither 
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the highly non-linear mapping of time onto the HR diagram nor the non-uniform mass and metallicity distributions of 
the stars in the galactic disk. Consequently, the derived age distribution is biased toward an older age compared to the 
real distribution. To avoid this bias, one needs to account for the a priori distribution functions of stellar parameters. 
For instance, the longer main-sequence timescale of lower-mass stars results in a smaller likelihood of observing low- 
mass post-main-sequence stars relative to higher mass stars. Bayesian probability theory including physically motivated 
prior distribution functi ons has been demonstrated to be an effect i ve means of determining unb iased ages of stars in 
the solar neighborhood (|Ng fc Bertellilll998t ILachaume et~aTlll999l: I"j0rgensen fc Lindegrenll2005t ). 

To model the SPOCS stars, we have constructed large and fine grids of theoretical stellar evolutionary tracks, 
computed with the Yale stellar evolution code (YREC). With YREC, the physical structure of the star can be calculated 
for each snapshot in time. Our procedure for evolutionary track calculations with YREC are described in detail in § 13.11 
Tracks have been computed for more than 250000 stars with slightly different initial conditions in mass, metallicity 
and helium abundances (see § 13.21 for the full description of the grids of stellar tracks) . The high resolution of the 
grids in time (lMyr) can apply accurate weighting for the accelerating evolutionary phases from the main sequence 
to post-main sequence. Each of the stellar models from each of the stellar evolutionary tracks is then assigned to a 
cell in a four dimensional grid of observable stellar parameters to increase the computational efficiency. Subsequently, 
posterior probability distribution functions (PDFs) are derived for each stellar parameter, in the framework of the 
Bayesian probability theory account (§ 13.30 . 

(iv) PDF error estimate — A calculated PDF is normally summarized by a single value for the best estimate 
and/or a credible interval. However, the choice of the best-fit value and the associated credible interval is often 
rather arbitrary. There are at least three choices of statistical values commonly used to summarize the PDF: the 
mean, median and mode. Since derived PDFs are not necessarily unimodal, let alone normal distributions, a single 
summary statistic will not always be sufficient to accurately describe the post erior PDF. Thorough analyses of accurate 
parameter representations have been done bv .Jorgens en fc Lindegrenl (|2005h with extensive Monte Carlo simulations. 
They have introduced the notion of "well-defined" age to distinguish the derived paramete rs whose posterio r PDF 
have a well-defined peak within the given parameter range. The same notation is adopted by I Nordstrom etaD (l200l 
to derive the ages for more than 10000 stars from the Geneva-Copenhagen survey. We hav e carefully presented all 
the ne cessary information from the calculated PDFs, applying a similar notation to that by Jorgcnse n fc Lindegrenl 
( 2005) . We will discuss the details of parameter estimates in § 13.4.11 

This paper is organized as follows. In § 2, we briefly describe the stellar sample, observations, data reduction and 
spectrum synthesis procedures. In § 3, we describe the stellar evolution code (§ 13. 1|) . construction of grids of stellar 
evolutionary tracks (§ 13.2ft , the calculation of derived stellar parameters in the Bayesian framework (§ 13.3ft and the 
error analysis for these parameters (§ 13. 4ft . The derived stellar properties and the parameter distributions for the 
SPOCS catalog stars are presented in § 4. The newly determined stellar properties of five planet-host stars previously 
modeled by FRS99 are presented in § 15.11 In § I5.2( we select five planet-host stars with particularly interesting stellar 
or planetary properties and discuss the derived models for these stars. The overall stellar parameter distributions 
are further analyzed in § I5.3( for parameter correlations and constraints on various dynamical formation scenarios for 
extrasolar planets. 

2. THE SPECTRAL ANALYSIS TECHNIQUE 

VF05 have carried out a uniform spectroscopic analysis of 1944 spectra for a sample of 1140 FGK stars in planet 
search programs at Keck Observatory, Lick Observatory and the AAT. The stars for these surveys are selected to 
optimize the achievable Doppler velocity precision and favor bright, chromospherically inactive, main-sequence or 
subgiant stars (My > 3.0, V < 8.5 and B — V > 0.5). Known stellar binaries with separations less than 2 arcseconds 
were rejected because the presence of a close stellar companion complicates the Doppler analysis. 

A detailed description of the methodology and an assess ment of random and system atic errors is provided in VF05. 
The spectral synthesis modeling program, SME (VF05, IValenti fc Piskunovl 119961 ) ass umes local th ermodynamic 
equilibrium and drives a radiative transfer code u sing Kurucz stellar atmosphere models (lKuruczHl993l) . and atomic 
line data (Vienna Atomic Line Database, VALD, iKupka et al.lll999t iRvabchikova et al"1ll999t ) to create a synthetic 
spectrum. The code employs a non-linear least squares Marquardt fitting algorithm to vary free parameters (T e g, logg, 
v sin i and abundances) in order to best match continuum and spectral line profiles in selected wavelength regions of an 
observed spectrum. With each Marquardt iteration, the SME program interpolates over a grid of 8000 Kurucz stellar 
atmosphere models before generating a new synthetic spectrum. Most of the SME analysis was made for wavelengths 
between 6000 — 6200 A to minimize problems with line blending. An additional wavelength segment 5175 — 5190 A 
was included to leverage the gravity sensitivity of the Mgl b triplet lines. 

As discussed in VF05, the ler uncertainties in the spectral modeling correspond to about ±44 K for T e g, ±0.06 dex 
for logg, ±0.03 dex for abundances, and ±0.5km/s for vsini. A number of researchers have carried out spectroscopic 
analyses for smaller subsets of stars in common with the VF05 sample. Effective temperatures and abundances show 
excellent agreement with these published results, however, small systematic offsets are seen. Comparisons of log g values 
show larger RMS scatter, but generally agree within quoted (~ 0.15 dex) uncertainties. While different spectroscopic 
analyses show reasonable agreement, the non- negligible systematic offsets demonstrate that spectroscopic results from 
different investigators should be combined with caution, particularly when looking for subtle correlations or trends. 
This large homogeneously analyzed sample is ideal for comparisons to the stellar evolutionary models here. 
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Fig. 1. — Sample stellar evolutionary tracks with the solar metallicity (X fa 0.71, Z fa 0.02) computed with YREC. Here 99 stars with 
known planetary companions (red solid circles) and 975 stars with no detected planetary companions (open circles) are overlaid. The 
thin solid lines with X are the theoretical evolutionary tracks of stars for every O.IMq from 0.5 to 2.0 Mq, evolved up to 14Gyr. The 
evolutionary tracks for stars with M < 1.0 Mq reach 14Gyr during the main-sequence phase. Each X on the tracks is separated by 20Myr. 
Note that the evolutionary sequence accelerates in the post-main-sequence evolution, corresponding to the Hertzsprung gap. The thick 
solid line represents the location of the ZAMS stars. The diagram represents only a small part of our grid of evolutionary tracks, which 
encompasses much wider ranges of T e g, L, [Fe/H], densely covered by the tracks that evolved from different initial conditions. 



3.1. Yale Stellar Evolution Code 

We use the Yale Rotational Evolution Code (YREC) in its non-rotating mode to calculate stellar m odels. YREC is 
a Henyey code which solves the equations of stellar structure in one dimension dGuenther et al.lll992f). The ch emical 
composition of each shell is updated separately using the nuclear reac tion rates of iGruzinov fc Bahcalll ([1998). The 
initial chemical mixture is the solar mixture of iGrevesse et all (|1993f ). scaled to match the metallicity of the star 
being modeled. Gravitational settling of helium and heavy elements is not include d in these models. For regions of 
the star which are hotter than logT(K) > 6, we use th e OPAL equat i on of state (|Rogers et al.l 1996V For regions 
where logT^) < 5.5, we use the equation of state from lSaumon et all (|1995l ). which calculates particle densities for 
hydrogen and helium including partial dissociation and ionization by both pressure and temperature. In the transition 
region between these two temperatures, both formulations are weighted with a ramp function and averaged. The 
equation of state includ es both radiation pressure and electron degeneracy pressure. We use the OPAL opacities 
(|Iglesias fc Rogers! [ 1996) for the interior of the star down to temperat ures of log T(K) — 4. For lower temperatures, 
we use the low-temperature opac ities oflAlexander fc Fergusonl (|1994[ ). For the surface boundary condition, we use 
the stellar atmosphere models of I Allard fc Hauschildt (1 19951). which include molecular effects. We use the standard 
Bohm-Vitense mixing length theory (|Cox fc Giulilll968t lBohm-Vitenselll958D with a=1.7161. This value of a, as well 
as the solar hydrogen abundance, X Q — 0.70785, is obtained by calibrating models against observations of the solar 
radius (6.9598 x 10 10 cm) and luminosity (3.85 1 5 x 10 33 erg/s) at the present age of the Sun (4.57 Gyr). 

3.2. Constructing the Grids of Stellar Evolutionary Tracks 

Our grids consist of 246000 theoretical stellar evolutionary tracks, separated by small intervals of mass and metallicity. 
In this section, we describe the iteration method we have adopted to construct these grids of model tracks. 

First, we have constructed models of pre-main-sequence (PMS) stars with masses ranging from 0.5 to 2.0 M@ for 
every 0.1 M & . These PMS stars are modeled by solving the Lane-Emden equation for a polytrope of index n = 1.5. 
Each star in the grids has evolved from a polytrope with the mass closest to the desired value. Since our grids have 
much finer mass bins (dM = 0.001M©) and a wide range of metallicity (from -1.0 to +0.60 dex in [Fe/H]), we first 
rescaled the selected polytrope to the desired values of mass and composition (X,Y and Z), then numerically relaxed 
it. After these parameters are rescaled, the star begins its PMS evolution through the zero-age main sequence (ZAMS) 
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TABLE 1 

Resolution of the stellar model grid 



X 


range 


AX 


Mass 


0.5 - 2.0 M Q 


0.001 M Q 


[Fe/H] 


-1.0 - 0.6 dex 


0.04 dex 


dY/dZ 


0.0, 1.5, 2.5, 3.5 




age 


0.0 - 14.0 Gyr 


20 Myr (M < l.OM ) 






1 Myr (M > 1.0M Q ) 



into the main sequence (MS) and subsequent evolutionary stage. The tracks through the PMS phase up to the ZAMS 
are not included in our modeling, since the SPOCS catalog would reject such young stars on the basis of strong Ca H 
and K emission 

Using YREC we calculate all the model parameters of interest (L,T e g, log g,R and so on) for each snapshot of 
the evolutionary sequence, starting from the ZAMS, until either (i) 14 Gyr is reached, or (ii) the end of the subgiant 
branch (beginning of the red-giant phase), whichever comes first. Table [1] describes the sets of initial stellar parameters 
from which stars are evolved. We vary the mass from 0.5 to 2.0 Mq for every 0.001 Mq, and [Fe/H] from -1.0 to 0.6 
dex for every 0.04 dex. YREC takes the compositional mass fractions {X, Y, Z} as input, thus the iron abundance 
is conve rted into the mass ra tio of heavier elements, Z, as Z = Zq x 10[ Fc / h 1, using the assum e d sola r value Zq — 
0.0188 (|Grevesse et al.|[T993D . For the helium abundance Y, the previous work b y iFord et al.l (fl999l) adopted the 
canonical linear relation betw een the abundances of helium and heavier elements (|Pagel et al.l fl992; Brcssar Tet al.l 
U994tlEdvardsson et al.lll993allrj ). 

- = (1) 

6Z~Z-Z Q 7 [) 

applying the solar helium abundance of Yq = 0.27335. We have also adopted the linear correlation between Y and Z, 
varying the values to be SY/SZ = 0.0, 1.5, 2.5 and 3.5, assuming a uniform prior distribution of 5Y/SZ. 

YREC solves a set of linearized equations using adaptive time steps St which are consecutively determined from a 
number of criteria. While YREC uses adaptive time steps to evolve the stellar models, we have interpolated all the 
evolutionary tracks to uniform time steps so as to properly account for the likelihood of observing a star at a given 
point in its evolution. Specifically, we have adopted At = 20 Myr for the stars with M < I.OMq and At = 1 Myr for 
the stars with I.OMq < M < 2.0M Q . The finer time resolution for the stars with higher mass is justified as follows. In 
the HR diagram, an evolutionary track of a star with a mass larger than ~ 1.2M Q has a sharp hook at the turnoff from 
the main-sequence track when the convective core contracts sharply due to the flat hydrogen profile at core hydrogen 
exhaustion (see Figure [1}. To account for this sharply non-monotonic behavior, a greater resolution of time is required 
for the stars with larger masses. When computing posterior PDFs, we weight each model by inverse of the time step, 
so as to account for the difference in the interpolated time step and to maintain the uniform prior distribution in age 
(see 



3.3. Bayesian Analysis 

To explore quantitatively the observational constraints on the ste l lar pa rameters, we employ the techniques of 
Bayesian inference. Our approach closely follows that of iPont fc Everl (|2004h . but we have generalized their methods 
as described below. Most notably, we: 1) include observational measurements of the surface gravity, 2) take into 
account the correlations between the derived values of the stellar metal abundance, surface gravity, and effective 
temperature, and 3) compute a dense set of stellar models that eliminates the need for interpolation between stellar 
tracks. 

In the Bayesian framework, the model parameters are treated as random variables which can be constrained by 
the actual observations. Therefore, to perform a Bayesian analysis, it is necessary to specify both the likelihood 
(the probability of making a certain observation given a particular set of model parameters) and the prior (the a 
priori probability distribution for the model parameters). Let us denote the model parameters by 9 and the actual 
observational data by d, so that the joint probability distribution for the observational data and the model parameters 
is given by 

p(d,6)=p(6)p(d\9)=p(d)p(6\d), (2) 

where we have expanded the joint probability distribution in two ways and both are expressed as the product of a 
conditional probability distribution (p(d\9) or p(9\d)) and a marginalized probability distribution (p(6) or p(d)). A 
marginalized probability distribution can be calculated by simply integrating the joint probability distribution over all 
but one of the variables (e.g., p{9) = J d9p(d, 9)). 

Often the model parameters contain one or more parameters of particular interest (e.g., the stellar mass, radius, and 
age) and other "nuisance parameters" that are necessary to adequately describe the observations but that we are not 
particularly interested in measuring for their own sake (e.g., the distance, helium abundance). 

We assume that the function mapping the stellar models to the observational data is a surjection or "onto function", 
meaning that we assume that our domain of stellar models includes at least one model that would result in any possible 
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set of true values of the observable stellar parameters. We use Bayes' theorem to calculate posterior probability 
distributions for the model parameters, as well as other physical quantities derived from the stellar models. 

Ideally, we would employ a hierarchical Bayesian analysis that could also incorporate the uncertainties in the the- 
oretical models (e.g., choice of equation of state, opacity tables, treatment of convection, etc.), so as to estimate the 
uncertainties in physical parameters even more accurately. Unfortunately, this is not yet computationally practical. 
While we have attempted to minimize such complications (at least for comparisons between stars in our sample) 
by making use of a single set of state-of-the art stellar models and the largest available uniform set of high quality 
spectroscopic observations and determinations of stellar atmospheric parameters, we acknowledge that our uncertainty 
estimates may not fully account for the systematic difference between different stellar evolutionary codes. 

3.3.1. Priors 

We assume a prior of the form p(M, r, D, X, Y, Z) = p(M)p(T)p(D)p(X , Y, Z). Thus, the priors for the mass, age, 
and distance are independent of each other as well as the chemical composition. However, the hydrogen, helium, 
and metal abundances are correlated. If we write p(X,Y,Z) = p(Z)p(Y\Z)p(X\Y, Z), then the logical constraint 
that X + Y + Z = 1, implies that p(X\Y,Z) = 5(1 — Y — Z), where 6 is the Dirac delta function. Given the 
observational difficulties in measuring Y for main-sequences stars, we base the helium abundance on the metallicity 
using an assumed value of 8Y/5Z. Thus, we assume p(Y\Z) = S(Y — F — (Z — Zq)SY/SZ). While chemical abundance 
studies of a variety of stellar populations suggest 6Y/8Z ~ 2.5 (see § I3.2[) . we employ a hierarchical model in which 
we assume a prior for SY/5Z ~ U[0, 3.5]. We write the prior for the metal fraction in terms of the surface metallicity, 

z = z Q io( [Fc/HI - [Fo/H1 °). 

Since the SPOCS catalog is intentionally biased towards high metallicities, we have chosen to adopt a uniform 
prior for the metallicity, p([Fe/H]) ~ U[— 1.0, 0.6], instead of using an empirical distribution for stars in the solar 
neighborhood. We implicitly assume that the abundance of all metals is proportional to the iron abundance. Thus, 
our prior can be written as p(M, r, [Fe/H] , D). 

For the prior for stellar mass, we use a truncated power law based on empirical estimates of the IMF, p(M) ~ M -2 35 
for 0.5Mq = M m i n < M < M max — 2.0M Q . Fortunately, the observations typically provide a tight enough constraint 
on the stellar mass that the results are relatively insensitive to the exact form of the mass prior. Since we impose 
sharp upper and lower cutoffs on the prior for the stellar mass, we are able to compute an extremely dense grid of 
stellar evolution tracks for relatively small range of stellar masses. The choice of the upper and lower limits is based 
on the selection criteria for the VF05 sample that we analyze in this paper. While our models and methods can be 
applied to other observations, the current set of stellar evolution tracks is limited to stars that are almost certainly 
between M m \ n and M max . 

We adopt a prior that is uniform in stellar age, p(r) = l/r max for yr < t < t max = 14Gyr. This choice represents 
the maximum entropy prior satisfying the obvious logical c riterion that all sta rs have a positive age that is less 
than the age of the universe from cosmological observations (jSpergel et al.l l20(33h . Since det ermining stellar ag es is 
notoriously difficult and the interpretation of the ages of a population of stars is subtle (e.g.. lPont fc Everll2004( l. we 
intentionally avoid incorporating prior observational or theoretical notions about the star formation history of the solar 
neighborhood. Indeed, we believe that this work has the potential to shed light on the history of star formation in the 
solar neighborhood, based on the combination of a large uniform sample of stellar atmospheric parameters determined 
with high resolution spectroscopy, a large, dense set of stellar evolutionary tracks, and our use of Bayesian inference. 
Note that this choice would be optimal for a stellar population that had a constant rate of star formation over r max . 
This is reassuring, since the observations are reasonably approximated by a globally constant star formation rate for 
the galactic disk. 

For the prior distribution of stellar distances, we assume a uniform density, p(D) ~ D 2 , for lpc = D m [ n < D < 
D max = lOkpc. The lower (upper) limit is intentionally chosen to be sufficiently small (large) that they are clearly 
excluded by the Hipparcos parallax measurements for all the stars that we analyze in this paper and hence do not affect 
the posterior probability distributions. The parallax is related to the distance by the definition 7r/arcsec = pc/D. 

Note that the above priors can be thought of as implicitly defining a set of priors for observables (e.g., parallax, 
luminosity, effective temperature, and surface gravity) via the mapping provided by the stellar evolutionary models 
(see § I3.1j) . We emphasize that a uniform age prior maps into a very non- uniform distribution in the observable 
quantities due to large changes in the time derivatives of the observable quantities with stellar age. Thus, we expect 
that our proper Bayesian treatment may result in significantly different and statistically superior estimates of physical 
quantities when compared to traditional frequentist analyses. Indeed, this is one of the primary motivations for our 
reanalysis of the VF05 sample. 

3.3.2. Likelihood 

We regard the observed parallax as being normally distributed about the actual parallax with a dispersion given 
by the standard error reported in the Hipparcos catalogue. In practice, the uncertainty in the stellar luminosity is 
dominated by the uncertainty in the parallax measurement, so the visual magnitude and bolometric correction can 
be regarded as exact. Note that the combination of parallax, visual magnitude, and bolometric correction place an 
observational constraint on the stellar luminosity that is asymmetric and has a positive skewness. 

We assume that the observational uncertainties in the stellar visual magnitude and parallax are independent of each 
other and the derived atmospheric parameters. However, we do account for the significant correlations between the 
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derived effective temperature, chemical composition, and surface gravity. We assume that the stellar atmospheric 
parameters derived from spectroscopic observations by VF05 are distributed according to a multi-variate normal 
distribution about their true values. Unfortunately, adopting a covariance matrix (C) based on the information matrix 
evaluated at the best-fit set of stellar parameters leads to significantly underestimating the observational errors, as 
determined by comparing analyses of the same star using multiple spectroscopic observations. This is due to the % 2 
surface being very "bumpy" and the curvature at any one local minima not reflecting the probability of the true solution 
being in another local minima. Since VF05 have multiple observations of several stars in their sample, they are able to 
compare the atmospheric parameters that they derive from multiple observations of the same star to more accurately 
determine their measurement errors (see VF05). We extend this approach to also determine the off-diagonal terms 
of the covariance matrix based on 9 sets of the derived solar atmospheric parameters (based on multiple observations 
of Vesta). For each star, we scale each of the empirically determined solar covariance terms by the standard error of 
both the relevant parameters for the target star. 
Thus, our likelihood function is 

L(V,tt, log T eff , [Fc/H] ,log. 9 )=i H (U,7r,logT off , [Fe/H] , log g) x L s (V,n,logT eS , [Fe/H] , log <?) 



—, — rrTo x ex P 

(27r) 3 / 2 a T 

1 

x — = cxp 

v /2 7 rdct(C) 



(7r obs - 10-( v ob 3 -Mv., mod )/5+n 



2 V^obs - ^mod) C (d obs - d mod ) 



(3) 



where a v is the measurement uncertainty in the parallax, dct(C) is the determinant of the covariance matrix adopted 
for <if bs , Mv,moA — Mv,mod(M, r, [Fe/H] , 5Y/SZ) is the absolute magnitude in the V band from the model (that 
includes a bolometric correction), and d^ od — d^ od (M, r, [Fe/H] , 5Y/5Z) is the set of logT e ff, [Fe/H], and logg for the 
model of a star of mass M, metallicity [Fe/H], age r, and helium abundance implied by 5Y/SZ. Note that we have 
separated the likelihood into two components, one that is a function of V a h s and 7r b s , and another that is a function 
of the spectroscopic parameters, d^ bs . 

3.3.3. Posterior 

The posterior probability for a set of model parameters (M, r, [Fe/H] , D) given the observational data d is given by 

p(M, r, [Fc/H] , SY/SZ, D\V ohs , n ohs , logT off , obs , [Fe/H] obs , log 3obs ) = 

p(M, r, [Fc/H] , SY/SZ, D)L(V ohs , ^ obs , r eff , obs , [Fe/H] obs , log 5obs ) 



/ dV dir dT eS d [Fe/H] d\oggp(M, t, [Fe/H] , 5Y/SZ, D)L(V, tt, T eff , [Fe/H] , log g) ' 



(4) 



where the integral is formally over the entire range of possible visual magnitudes, parallaxes, effective temperatures, 
metallicities, and surface gravities. When there are multiple spectroscopic observations of the same star, then the 
spectroscopic portion of the likelihood, Ls, is replaced by a product of multiple Ls's, with one Ls being evaluated 
with each of the observed sets of spectroscopic parameters. 

Often, we are particularly interested in the posterior probability density marginalized over all but one of the model 
parameters. This can be easily calculated from the posterior probability density by integrating over all but one of the 
model parameters. For example, the marginal density for the stellar mass is given by 

p(M\V ohs , 7r obs , logT off!0bs , [Fe/H] obs , log 5obs ) = 

dr d [Fe/H] dSY/SZdDp(M, r, [Fc/H] , SY/SZ, D\V ohs , 7r obs , logT cff , obs , [Fe/H] obs , log. 9obs ). (5) 



/ 



We directly calculate the marginal posterior densities for the stellar mass, age, and metallicity. We also calculate 
marginalized posterior densities for derived physical quantities. For example, the marginal posterior density for the 
stellar radius, R, is given by by 

p(R\V ohs , 7r obs , logr of f !0bs , [Fc/H] obs , log.g obs ) = 

dM dr d [Fe/H] d5Y/SZ dD S(R - R mod (M, r, [Fe/H] , SY/SZ)) 

p(M, r, [Fe/H] , SY/SZ, D\V ohs , 7r obs , logT off , obs , [Fc/H] obs , log 3obs ), (6) 

where R mo d(M, t, [Fe/H] , SY/SZ) is the radius of the stellar model with mass M, age r, metallicity [Fe/H], and helium 
abundance implied by SY/SZ. 

3.3.4. Numerical Methods 

The main difficulty in performing Bayesian analyses is the difficulty of performing all the necessary integrals, par- 
ticularly in high dimensional parameters spaces. Here, we describe the numerical methods we use to approximate the 
necessary integrals. 
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To numerically compute these marginal densities, we discretize the integrals in the standard way. So Eqn. [6] becomes 
p(R <R<R a + AR\V ohs , 7r obs , logT cff , obs , [Fe/H] obs , log 5obs ) S£ 



i j k I 

I ' 



Ro < Rmod [ M h T fe , [Fe/H^ , ^- I < R Q + AR 



5Zi 
SY 

xp(Mi,T k , [Fe/B]j , — ,£>|Vob8,7rob s ,logT eff , ob8 , [Fe/H] obs , log^obs), (7) 

where AM, is the spacing between the ith and i + 1th stellar mass in our grid of stellar models, A [Fe/H]^ is the spacing 
between the jth and j + 1th metallicity in our grid of stellar models, A|^ ; is the spacing between the Ith and I + 1th 

value of SY/5Z in our grid of stellar models, and Arfc is the spacing between the fcth and k + 1th age is the spacing of 
the outputs in the stellar track with mass Mj and metallicity [Fe/H]^. The indicator function I(R < R < R + AR) 
is defined to be 1 if R < R < R + AR and otherwise. For values of the various spacings between parameter values 
in our grid of stellar models, see Table 1. The stellar evolutionary tracks are computed with a variable time step, so 
as to efficiently evolve the star rapidly during the main sequence, but provide the necessary temporal resolution to 
accurately model the early and late stages of evolution. In order to facilitate the numerical integration, we interpolate 
within each stellar track to provide each of the observable and derived physical quantities at a series of stellar ages, 
each separated by Ar^ = At. The value of At varies between stars to reflect the speed of stellar evolution for each 
track (see Table 1). It is important to note that we interpolate only in time within a single evolutionary track, and 
not between stellar tracks of different masses or compositions. Thus, our interpolation does not suffer from any of 
the complications typically associated with stellar isochrone fitting. For a given set of observations, we can compute 
the marginal quantities by summing the product of the prior times the likelihood evaluated at each time of each 
evolutionary track. 

In practice, the measurement uncertainties for each of the parameters is much smaller than the entire allowed range of 
the parameter. Therefore, we can approximate each of the above summations over all model parameters by summations 
over a region R' that contains all the points that contribute significantly to the integrals. We are quite conservative, 
and choose R' to include every point for which |log T e g , bs - logT offimod | < 10er Tcff , |log [Fe/H] obs - log [Fe/H] mod | < 
10a [Fc/H] , | log g bs - log swi | < 10<r log5 , and |^ obs - lO"(^- M v, m o d (M.r,[Fc/H]))/5+i | < ^ We manually check 
that reducing the volume of R! results in no significant difference for the marginal distributions. 

3.4. Characterizing the Derived PDFs 

While we will provide the full posterior distributions upon request, it is often convenient to summarize the posterior 
PDFs. Calculated PDFs are typically represented by two quantities — a single "best estimate" value, and some 
associated credible interval. However, the derived PDFs can often manifest complicated shapes that are not readily 
fit by a standard normal distribution (for example, see Figure [2] for sample age-PDFs). Thus, defining the best 
representative value and determining the credible interval is a non-trivial task. An incorrect procedure for model 
parameter estimation can introduce an unwanted bias in a large sample and also fail to extract all the necessary 
information from the PDFs. A use f ul rep resentation of derived age probability distributions has been thoroughly 
discussed bv lJorgensen fc Lindegrenl (|2005h . 

3.4.1. Best Estimate 

The main goal of this section is to define a single estimate value x, along with the range of plausible values [xi, X2] 
determined from the selected credible level. Commonly used statistical quantities for representing the best estimate 
from the PDFs are the median (the bisector of the area under the PDF curve), the mean (the expectation value), and 
the mode (the most probable value). No matter which statistical quantity is selected, it inevitably includes certain 
arbitrariness and statistical biases. Thus, when studying the statistical properties of a large sample of stars, one needs 
to be fully aware of the type of biases introduced and separate them carefully from the true values of interest. 

We have chosen the mode, i.e., the global maximum of the derived PDF, as our primary estimate of the model stellar 
parameter (denoted by the hat — e.g., f for the best age-estimate). The mean and the median have their advantages 
in that they always reside within any given parameter range. For instance, the maximum age of HD10700 in Figure[2] 
apparently lies beyond 14Gyr, yet it is still possible to present the estimated age as 11.6 Gyr by choosing the median 
value. However, as pointed out by Jorgcnsen & Lindegren (2005), when the derived PDF is broad due to a poor 
observation, the mean or the median tends to be deviously located near the middle of an arbitrary parameter range. 
For example, in Figure [2] it is in no sense meaningful to say that the estimated age of HD144628 to be 7.5 Gyr (the 
mean), or 7.7 Gyr (the median). It is evident that these estimators would lead to a spurious bias toward 7 Gyr, the 
central value of our selected age interval, 0-14 Gyr, when analyzing a population of similar stars. Furthermore, the 
calculated mass or age-PDFs often have multiple peaks as seen for HD193307 in Figure[2j since the model parameters 
(M,t, Z,etc) are not in one-to-one relation with the observed parameters (My, T e g, [Fe/H]). Multiple observations 
for a single star can lead to bimodal posterior distributions that are even more difficult to summarize with a single 
estimator. 
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Fig. 2. — Normalized age probability distribution functions (PDFs) of sample stars. HD101614 (solid curve) and HD193307 (dot-dashed 
curve) have "well-defined" ages. HD193307 has a secondary maximum in the age-PDF, T2 = 3.4 Gyr with a relative probability of 23% 
with respect to the best-estimate age, f = 6.5 Gyr. Only the best-estimate age and the upper-bound of the 68 % credible interval for the age 
of HD120237 (short-dashed curve) will be specified in the table (Thigh = 2.9 Gyr). Only the lower-bound for the age of HD10700 (dotted 
curve) will be specified (ti ow = 12.2 Gyr). No meaningful age can be derived for HD144628 (long-dashed curve) . The credible intervals are 
defined as the range between the points where the PDFs cross the solid horizontal line, P(r) = 0.6 (see S I3.4,2( . 

3.4.2. Error Estimate 

The conventional definition for ±lcr credible interval is to assume a Gaussian distribution and find the 
range [iBi,^] such that the fractional area under the curve falls to 68% of the total integrated probability; 

P{x)dx/ JJ max P(x)dx — 0.68. Instead, we normalize the PDFs such that P(x) = 1.0 and estimate our "la- 
error" [x\, X2] to be the interval between the two points where P(%i) = P(x2) — 0.6 (jNordstrom et al1l2004[ ). This is 
a more appropriate choice for age-PDFs, since for many stars age is a particularly weakly constrained parameter. For 
example, the derived age-PDF for HD193307 in Figure[2]has a smaller secondary maximum. For such ambiguous cases, 
the credible interval estimated by the fractional area under the curve does not describe the correct uncertainty range 
for the best age estimate (the credible interval may even lie outside the mode value). Also, computing the la credible 
interval from 68% fractional area is likely to underestimate the uncertainty of broad age-PDFs (e.g., HD144628). Note 
that for a standard Gaussian distribution, the credible interval [xi,^] defined as P(xi) = P[x2) = 0.6 is the region 
in which the true parameter lies for 68% of all the cases. This defin i tion o f a credible interval can be also applied to 
non-Gaussian type model parameter PDFs. iJorgensen fc Lindegrenl (|2005| ) have done extensive Monte Carlo simula- 
tions using 10 3 synthesized stars with typical observational errors, and confirmed that 68% of the recovered ages fell 
within [21,2:2] of the true age. 

3.4.3. Our Notation for the Parameter Estimates and the Credible Intervals 

Here we summarize the notations we have used to characterize the parameter estimates and credible intervals in 
Table [2] When the credible interval around the m ode is entirely within th e parameter range of our grids, [a; m ; n , £ max ] , 
the estimate is called a "well-defined" parameter (jNordstrom et al1l2004f ). However, it often happens, particularly in 
the age-PDFs, that the credible interval is truncated on one or both sides of [x m i n , x max ] (e.g., HD120237, HD10700 
and HD144628 in Figure^. In these cases, the truncated side of the credible interval is left blank in the table. This 
means that only an upper or lower bound (or, possibly, no bound at all) can be specified at the 68% credible level for 
the considered parameter. Similarly, sometimes the mode is not well-defined, as the position of the global maximum x 
coincides with either boundary of the permitted parameter range (e.g., P{x raax ) is unity for HD144628). In this case, 
the best-estimate value i is not specified in the table. However, an upper or lower bound of the credible interval may 
still be specified, if it exists. 

Lastly, the calculated model parameter PDFs often exhibit a bimodal feature (e.g., HD193307). Sometimes the 
secondary maximum lies outside the credible interval of the primary. In these cases, the secondary maximum and its 
probability relative to the global maximum are also specified in the table. 

It is often useful to also provide summary statistics for each of the model parameters. By their very nature, summary 
statistics throw away much of the information that is contained in the posterior PDFs. The loss of information is 
particularly acute for complex PDFs (e.g., very broad, highly skewed, multi-modal), such as often occur with the 
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TABLE 2 
Theoretical Stellar Parameters 



Stellar ID 




M 


Mow 




M 2 




7~low 


Thigh 


T2 




R 






Rcc 




log 9 




[Mq] 


l M @\ 


l M &\ 


[Mq\(/o) 


[Gyr] 






[Gyr] 


(%) 




1-n-eJ 










[ C 9 S ] 


HD102158 





.914 


0.890 


0.936 




11 


.36 


9.84 


12.92 






1 


i i -0.05 
•• L1 +0.05 





n1 7-0.003 
■ ul ' +0.003 


n 01,-0.023 

U.JlO +0 023 


4 


q O -0.04 

""+0.04 


HD102357 


1 


.129 


1.105 


1.153 




3 


.48 


2.56 


4.20 






1 


qo-0.05 
■ ZJ +0.06 





nn o-0.001 
■ UUJ +0.001 


0.218 +0017 


4 


qo-0.04 
••"+0.04 


HD102365 





.889 


0.858 


0.923 




9 


.48 


6.44 


12.44 









Qfi-003 
■ yD +0.03 





n01 —0.003 
■ uzl +0.003 


(\ ryrjri — 0.018 

u - z ' z +0.018 


4 


4 ,-0.04 
'^°+0.03 


HD102438 





.868 


0.840 


0.903 




10 


.04 


6.40 


13.40 









QO -0.03 
■ yz +0.03 





026-0 ° 03 
■ UZD +0.003 


270-°' 019 
u - z,u +0.016 


4 


-7-0.03 
■^'+0.04 


HD102634 


1 


.322 


1.296 


1.356 




2 


.80 


2.60 


3.04 






1 


pi -"■07 





nn4 -o.ooi 
■ UU4 +o.ooi 


O 9Q9 — 021 
u - ZMZ +0.035 


1 


17-0.04 
' i '+0.03 


HD 102870 


1 


.353 


1.319 


1.381 


1.236 (0.25) 


2 


.96 


2.64 


3.20 


4.36 (0.27) 


1 


7n -0.05 
■ ,u +0.04 





006-0 ° 03 
■ UUD +0.002 


340-0' 031 
U.J4U +Q 022 


4 


14-0-04 
• i ^+0.04 


HD 102902 


1 


.739 


1.707 


1.883 




5 


.84 


5.32 


6.68 






2 


,,,-0.25 
■^°+0.04 





i r,n-0.033 

■ lz,;> +0.019 


u - ,UO +0.093 


3 


q ,-0.02 
' MO +0.05 


HD103095 





.661 


0.655 


0.689 








2.44 











■ DO +0.01 





n r, o -0.001 
■ uzz +0.002 


1 70-0 006 


4 


fi o-0.01 
• DJ +0.03 


HD103432 





.948 


0.918 


0.970 










3.64 









oq-0.03 
■ o9 +0.03 





n2fi -0.003 
■ UZD +0.0O3 


242-0 00 9 

u - z ^ z +0.013 


1 


,q-0.03 

' OJ +0.03 


HD10360 





.750 


0.732 


0.767 










0.60 









71 -0.03 
■ (1 +0.02 





n4f) -0.002 
■ u4u +0.004 


21 1 -0.013 

u - zll +o.on 


4 


.64-0-0 2 


HD10361 





.761 


0.746 


0.783 










0.52 









7o -0.04 
■ ,z +0.03 





n q7-0.001 
■ u,:s '+0.005 


205-0' 008 
u - zuo +0.013 


1 


64-0.02 


HD103829 





.984 


0.962 


1.006 




9 


.00 


7.28 


10.48 






1 


ln -0.07 
■ iu +0.08 





n24 -0.004 
■ uz *+0.003 


n oncr-0.030 
U - JZO +0.029 


4 


0,-0.05 

• ,3o +0.06 


HD103932 





.775 


0.769 


0.783 










0.68 









7 O -0.02 
■ ,z +0.02 





nfil -0.003 
■ UDi +0.003 


n -0.005 
u.zo± +0 006 


1 


fi q-0.02 

■ DJ +o.oi 


HD104067 





.831 


0.799 


0.845 







.40 




6.48 









■77-0.02 
■ ' '+0.02 





n , 2 -0.004 
■ uoz +0.002 


241 -0' 007 

u - Z41 +o.oio 


4 


fin -0.03 
DU +0.02 


HD 104304 





.980 


0.936 


1.024 




8 


.48 


5.68 


11.04 






1 


oi -°- 03 

■ ul +0.03 





oo, -0.004 
■ u ' 3o +0.004 


Q14-O.020 


4 


<q-0.03 

'^ J +O.03 


HD10436 





.624 


0.616 


0.628 










0.44 









■ 09 +o.oi 





n q 7 -0.002 
■ uo1 +0.006 


174-OOO2 
u - i '^+0.004 


4 


7n -o.oi 

• ,u +0.02 


HD104556 






1.897 




1.403 (0.36) 


2 


.72 


2.48 


3.36 






3 


-i o -0.05 
■ iz +0.05 





QH1 -0.048 
■ t " 5i +0.032 


1 , 7 n-0.031 
±.orz +0 lig 


3 


7fi -o.oi 

■ ,D +0.01 


HD 104576 


1 


.035 


1.011 


1.051 










0.52 









OK-0.03 
■ yo +0.03 





nlfi -0.002 
■ U±D +0.003 


n O qq-0.003 
u.zoo +0 006 


4 


,1-0.01 

• ol +0.02 


HD 10476 





.816 


0.805 


0.838 








8.84 











oq-0.02 
■ oz +0.02 





n4 , -0.003 
■ u *°+0.002 


0.266+8;gig 


1 


, 4 -o.oi 

' oyi +O.03 


HD105 


1 


.129 


1.113 


1.145 










0.60 






1 


nfi -0.03 
■ UD +0.03 





nnfi -0.002 
■ uuo +0.002 


n nntr-0.008 
u - zzo +0.006 


4 


4 ,-o.oi 

• 4o +0.02 


HD105113 




1.201 


1.159 


1.227 






5.04 


4.68 


5.72 








1 QO-O- 10 
J " aa +0.13 




007 -0 - 002 
u ' uu '+0.003 


0-482;°;°i 




qo-0.04 

■J-^+o.os 


HD105113B 




0.951 


0.937 


0.965 










0.60 








u -°°+0.02 




n n o 7 -0.002 


243-0 00 4 

u - z *' 3 +0.006 




4 tro-0.02 
4 -O»+0 .01 


HD105328 




1.210 


1.133 


1.238 






4.36 


4.04 


4.72 


5.84 


(0.46) 




1 46-0' 05 
i - 4D +0.05 




01 o-o- 002 

u - ulz +0.002 






4 21 -° 02 
^• zi +0.03 


HD 105405 




1.074 


1.042 


1.100 






5.32 


4.60 


6.08 








1 29-° 07 
1 - ZM +0.06 




004-0' 001 


250-0' 019 
u - zou +0.020 




4 27-0- 05 
^• z '+0.03 


HD105631 




0.947 


0.921 


0.969 










2.80 








n o 7 -0.03 
U - 8 '+0.03 




041 -° 002 
u - u4i +0.003 


n 9 , 7 -0.008 
u - zo '+0.011 




4 56-° 03 
4 - oo +0.02 


HD106116 




0.981 


0.955 


1.013 






8.80 


6.76 


10.56 








1 06-° 04 




029-0' 003 
u.uzy +0 004 


n qr,r,-0.020 

u - c>zz +0.022 




4 QO-0.04 
4 - r>y +0.04 


HD106156 




0.957 


0.931 


0.983 










2.60 








- 88 +0!03 




042-° 003 
u - u4z +0.003 


260-o o° 9 

u - zou +0.011 




4 56-0' 03 
4 - OD +0.02 


HD106252 




1.007 


0.985 


1.031 






6.76 


5.16 


8.16 








, ,,-0.05 
1 - 11 +0.05 




o 01 5-0 00 3 

u - ulo +0.003 


O oqq— 021 
u - zao +0.020 




4 36-0- 04 
4 - JD +0.04 


HD 106423 




1.197 


1.141 


1.289 






3.68 


3.20 


4.12 








, sl -o.io 

1 -' 31 +0.10 




flu -0 ' 003 
u.ui4 +0 003 


r\ 001 — 0.025 
U - OZ1 +0.027 




4 - JU +0.05 



Note. — The entire catalog of derived stellar parameters is available through the California & Carnegie Planet Search web page, 
http : / /exoplanets . org/SPOCS^evol .html. 

marginal posterior PDF for the stellar age. In order to provide a succinct summary of the size of the tails of the 
posterior PDFs, we provide the median value as well as 68 %, 95 % and 99 % credible intervals for each stellar parameter 
in the electronic version of the table. 

4. DERIVED STELLAR PROPERTIES 

With the methods described in the previous chapter, we have computed the stellar parameter PDFs for 1074 SPOCS 
stars with T c g within the range 3700 - 6900 K. The derived stellar parameters for a small subset of the sample are 
presented in Table [2] The entire table of derived parameters is available in the electronic version of the paper and also 
in the California & Carnegie Planet Search website 6 . 

4.1. Stellar Ages 

Figure [3] shows the derived age distributions of stars in different sub-samples of the SPOCS catalog. The curves 
show the sum of the normalized age-PDFs of all the stars. Among the total of 1074 SPOCS stars we have used, 669 
stars have "well-defined" ages as defined in § 13.4.31 Among those, 606 stars have estimated errors smaller than 5 Gyr, 
and 169 stars have errors smaller than 1 Gyr. The large fraction of stars with the youngest ages (r < 1 Gyr) and the 
oldest ages (r > 13 Gyr) in Figure [3] are clearly artifacts from choosing the mode value as the best-estimate age. As 
discussed in § 13.4. 1[ the mode of a poorly constrained age-PDF tends to reside near either end of the selected age range 
(0 or 14 Gyr in this case). These accumulations at extreme ages are removed in the distribution of "well-defined" ages. 

Note that the selection criterion of the SPOCS catalog is such that precision of radial-velocity observations is 
optimized. The sample stars are selected by visual magnitude, but not the stellar distance. Conseq uently, the catalog 
contains more distant, early- type stars that are intrinsically more luminous. iFischer fc ValentH (|2005f ) defined a volume- 
limited sample with a radius of 18 pc, within which the number density of F-, G- and K-type stars per volume is nearly 
constant as a function of distance. Among the 1074 stars we selected from the SPOCS catalog, 203 stars are in the 
volume-limited sample. Figure [4] illustrates the bias toward intrinsically luminous stars in the whole sample. The 

6 http://exoplanets.org/SPOCS _evol. html 
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Fig. 3. — Distributions of derived ages. The right panel includes only the 99 stars in the SPOCS catalog with known planetary companions. 
The best-estimate ages are presented in the histograms, and the integrated total age probability distribution functions are presented as the 
curves. The blue dotted line includes only the stars with well-defined ages. The red dashed curve is the age distribution of the 61 stars in 
the volume-limited sample with well-defined ages. 




Fig. 4. — Distributions of the calculated luminosity for the whole sample (solid) and for the 203 stars in the volume-limited sample 
(dotted, blue). The volume-limited sample has a luminosity distribution peaked around logL/Z/Q = —0.32 whereas the luminosity of all 
the sample stars are peaked at log L/Lq = 0.15. A larger sampling volume includes more intrinsically bright stars. 



volume-limited sample is in fact rather abundant in intrinsically faint stars, whereas the whole sample consists of 
more stars with luminosity above the solar value. The age distribution of the volume-limited sample is slightly shifted 
toward a younger age, due to a larger fraction of relatively unevolved, less bright stars. 

The age distribution of stars with known planetary companions is presented in the right panel of Figure [3] The 
SPOCS catalog contains 99 planet-host stars, 74 of which have well-defined ages. The integrated age-PDF of all the 
stars has sharp local maxima, arising from the small number statistics and selection effects. The mode value of 3.3 Gyr 
of the planet-host stars coincides with that of all the well-defined ages, however, the peak is more distinct for the 
planet-host stars. The two maxima around 3.3 Gyr and 8.0 Gyr mostly consist of G-type stars, since planet-search 
programs typically target stars with solar-type spectra. Most of the stars in the 2 — 5 Gyr bins consist of young 
main-sequence stars with spectral type F8 - G4 and masses 1.0 - 1.9 M©, whereas the stars in the 7 - 9 Gyr bins are 
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Fig. 5. — Distributions of the derived masses. The right panel includes only the 99 stars in the SPOCS catalog with known planetary 
companions. The best-estimate masses are presented in the histograms, and the integrated total mass probability distribution functions 
are presented as the curves. The mass distribution for the stars in the volume-limited sample (< 18 pc) is presented as the blue dotted 
curves . 
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Fig. 6. — Distribution of the derived masses and the distances determined from the Hipparcos parallaxes. The red crosses represent the 
stars with known planetary companions. 

less massive G3 - G5 stars with masses 0.8 - 1.1M Q , many of them likely to be in the subgiant phase. 

Note that the stellar age is generally the most poorly-determined parameter. Nearly half of the SPOCS stars have 
derived age uncertainty greater than 5Gyr, which is comparable to the entire parameter range (0 - 14Gyr). Also, 
the age uncertainty is more sensitive to the stellar mass than on the accuracy of the spectroscopic observations (the 
derived age - mass relation is discussed in more details in ? I5.3.1[) . Figure [3] merely provides a graphical representation 
of the summary statistics and the contributions from different types of stars to the overall age distribution, and it 
should not be confused with the actual star-formation rate in the solar neighborhood. 



4.2. Stellar Masses 



Figure [5] shows the distributions of derived masses. The mass-PDFs are generally much better constrained than the 
age-PDFs. Well-defined masses are obtained for the majority (97%) of sample stars, including the ones with poorly 
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TABLE 3 
Properties of Sample Stars 



Observed Data 




Posterior 


Star T cff [K] [Fe/H] logg [cgs] 


T cff [K] 


[Fc/H] logg [cgs] 



t Boo 


6387±44 


+0.25±0.03 


4.25±0.06 


63907™ 


+0.31±0.04 


4.27 


51 Peg 


5787±44 


+0.15±0.03 


4.45±0.06 


5814-53 


+0.22±0.03 


4.36 


v And 


6213±44 


+0.12±0.03 


4.25±0.06 


61597. 42 


+0.16±0.04 


4.17 


55 Cnc 


5253±44 


+0.31±0.03 


4.45±0.06 


5327±49 


+0.37±0.04 


4.48 


p CrB 


5823±44 


-0.14±0.03 


4.36±0.06 


5855 + gi 


-0.17±0.05 


4.21 



-0.02 
+0.04 
-0.02 
+0.04 
-0.03 
+0.02 
-0.01 
+0.05 
-0.02 
+0.05 



Note. — The surface properties of the sample stars 
scopic values by VF05 and the posterior estimates. 



comparison between the spectro- 



TABLE 4 
Comparisons with Ford et al. 



Star 



Af» [M 



©J 



Age [Gyr] fl„[R, 



0J 



M ce [M ] Rco[R ] a logg [cgs] 



t Boo (HD120136) 
51 Peg (HD217014) 
v And (HD9826) 
55 Cnc (HD75732) 
p CrB (HD143761) 



1.34] 



-0.04 
"*+0.05 

1.37±0.08 
1.05±0.04 
1.05: 



1.64 



+ 0.44 



1.2 



-0.08 
J +0.09 
, 0,-0.01 
L - OL +0.02 



1.34] 



-0.12 
*+0.07 

u - ao +0.05 



6.76 
7.6 

3.12 
2.6 



+ 1.2 
1.48 



0.95 



-0.09 
+ 0.11 



0.96±0.02 



+ 1.64 
5.1 

+4.0 
0.24 
+ 0.20 
1.0 
+2.1 
> 7.24 

8-4+t? 



1.46±0.05 
i 41-0 09 

1.15±0.04 
1.16±0.07 

-0.05 
+0.04 
0.10 
+0.11 
0.01 
+0.03 
0.03 
+0.02 
0.02 
+0.03 
0.08 
+0.09 



1.64 
1.56 
0.93 
0.93 
1.35 
1.35 



002-°' 002 
u ' uuz +0.003 

< 0.002 

023-°' 004 

n n23-°- 006 
u.uzo +0 007 

005- 003 
u.uuo +0 002 

002- 002 
u.uuz +0 003 

u ' uo '+0.005 
n 046- 004 

015- 002 
u ' ulo +0.003 
n 0,00-0.009 



0.23±0.01 
n 99— 0.18 

u - zz +0.19 
u -' 3o +0.03 

0.34±0.11 

-0.07 
+0.03 
0.10 
+0.17 
0.03 
+0.01 
0.06 
+0.05 

0.38±0.02 
0-477_°,i 3 



4.27 
4.27 
4.36 



0.32 
0.27 
0.32 
0.30 



-0.03 
+0.04 
0.07 
+0.05 
0.03 
J +0.04 

4.33±0.09 

-0.04 
+0.02 
0.10 
+0.07 
-0.01 
-0.05 
-0.07 
+0.04 
0.01 
+0.03 
0.06 



4.16 
4.18 
4.48 
4.50 
4.18 
4.13 



Note. — Comparisons with the calculations from Ford ct al. (1999). Their results arc shown 
second row of each star. 

a Width of the convcctivc zone, measured from the outermost stellar surface 



the 



determined ages. The derived masses have a median of 1.03M©, consistent with the median of 1.07M Q quoted by VF05 
within the la uncertainty. Since only the stars with larger masses and luminosities are selected at large distances, there 
is a clear correlation between the derived mass and the stellar distance (see Figure • The volume-limited sample 
contains more stars with sub-solar masses, with a median value of O.87M0. 

The small gap in the distribution around 0.9M Q is likely to originate from our choice of age cutoff at 14 Gyr. The 
mass of 0.9M Q roughly corresponds the turn-off mass for this age (the critical mass beyond which a star can evolve away 
from the main sequence within 14 Gyr; see, for example, the sample evolutionary tracks for 55 Cancri in Figure ITTj). 
Indeed, the 125 stars in the mass range 0.8 - 0.9Mq have particularly old ages: 39 stars (31%) with ages > 12 Gyr and 
58 stars (46%) with > 10 Gyr. Note that in the complete sample only 7% of the stars have derived ages r > 12 Gyr, 
and 16% have r > 10 Gyr. 

The mass distribution of the 99 known planet-host stars is also presented in the right panel of Figure [5] The 
mass distribution of planet-host stars behaves similarly to that of the entire sample, with a median mass of 1.07M Q . 
The lower-mass cutoff near ~ 0.8M Q comes from the sampling criterion (V < 8.5) of the SPOCS catalog. Also, 
stars with masses below ~ 0.8M Q correspond to M-dwarfs, which are typically excluded from planet-search programs 
because of their low luminosities and strong chromospheric activities, which affect the radial-velocity measurements 
(|Delfosse et al.ll2~00Q[ ). All the planet-host stars with derived masses < 0.9Af© in our sample are KO - 3 stars. There 
are only four stars in the extrasolar planet catalog with spectra later than K3: HD63454 (K4V), GJ436 (M2.5), G1581 
(M3) and GJ876 (M4V). None of these stars is included in the SPOCS catalog. 

5. DISCUSSION 

5.1. Comparisons with Paper I. 

FRS99 modeled five star s with extrasolar planets k nown at the time, using the observed parameters obtained by 
lGonzale3 (|1997l 119981 ) and iGonzalez &: Vanturd (|1998l ) . In this section, we compare those results to our new models 
using the spectroscopic observations by VF05. The overall comparisons are summarized in TableJH and the theoretically 
derived atmospheric properties are listed in Table [H The normalized PDFs of the stellar properties are presented in 
Figure \7\ 

5.1.1. t Bootis (HD120136) 

A planet with a mass M p isini w 3.9Mj around the young F7V star rBoo was discovered bv lButler et al.l (I1997T) i n 
a 3.31-day period. Despite its assumed young age, the star shows little photometric variability ( Baliunas et al.lll997h . 
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Fig. 7. — The derived PDFs for mass (M), age (t), radius (R), mass and width of the convective zone (M ce , R cc ) and surface gravity (logg) 
of five planet- host stars: r Boo (solid, blue), 51 Peg (dark-green, short-dashed), v And (purple, dot-dashed), 55 Cnc (green, long-dashed) 
and pCrB (red, dotted). The derived parameters and credible intervals are listed in Table l4l 

Among the five planets discussed in this section, it is the hottest and the second most metal-rich star after 55 Cnc. 
The mean effective temperature T e g = 6410 K and metallicity [Fe/H] = +0.26 are derived by VF05 using the SME 
with Lick/Hamilton spec tra at three different epochs. T he planetary system rBoo is associated with a visual binary 
companion at ~ 240 AU (|Hald lT994; Pati ence et al.|[2002T ). The stellar companion does not affect the spectral analysis 
because of its large orbital separation. 

Our calculated properties of rBoo are consistent with the model made by FRS99: tBoo is a young, fairly mas- 
sive main-sequence star with a well-constrained age esti mate, 1. 64^44 G yr. The activity - age relation determined 
fro m the mean Ca I I flux suggests an age of ~ 2 Gyr (Bali unas et al.l 119971 ) . Th e rev ised spectrosco p ic an alysis 
by I Gonzalez fc Laws! (|2000l ) has been applied to isochrones bv ISchaller et alT(|1992f ) and ISchaerer et al.l (|1993l ) and 
yielded an age 3.3 ± 0.5 G yr. All the other isochrone or surface activity analyses estimate the age to be < 3.1 Gyr , 
consistent with our results (jFuhrmann et alj|1998t iLachaume et~aHll99fllSuchkov fc Schultzll2001l : [Henry et alj|2000at 
iNordstrom et al.ll2004h . 
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Fig. 8. — Theoretical HR diagram and stellar evolutionary tracks for r Bootis. The observed T c g from the SPOCS catalog and luminosity 
calculated from the observed magnitude and Hipparcos parallax are shown with la uncertainty. All the evolutionary tracks inside the 
dashed rectangle (lOcr observational uncertainty) are used for the stellar model calculation. The tracks are presented from ZAMS onward, 
with our best-fit metallicity [Fe/H] = +0.32. Each X on the tracks is separated by 100 Myr. Note that this time step is exaggerated in this 
figure to illustrate the acceleration of the evolutionary sequences — in the actual grids used for the calculations, the time resolution is as 
fine as 1 Myr. Theoretically derived model parameters of tBoo are M = 1.34Mq and age r = 1.64 Gyr. 



Figure [8] shows the theoretical evolutionary tra cks with [Fe/H] = +0.31 . This metallicity est imate is consistent wit h 
other spec troscopic observatio ns: +0.32 ± 0.06 ((Gonzalez fc Laws! 120001 ) and +0.27 ± 0.08 (jFuhrmann et al.l fl998). 
VF05 and lSantos et all (|2004f ) have derived slightly lower values, 0.25 ± 0.03 and 0.23 ± 0.07 , respectively. 

5.1.2. 51 Pegasi (HD217014) 

The first extrasolar planet around a solar- type star was discovered bv lMavor fc Que loz (1995), around the G5V star 
51 Pegasi. The planet has a mass M p isini w 0.47Mj and an orbital period of 4.23 days. 
From the three Lick/Hamilton spectra, VF05 derived T off = 5787 ± 44 K and [Fe/H] = 0.154 ± . 029, consistent 



with other PTE analyses, T fiff = 5750 - 5820 K and [Fe/H 



0.14 - 0.21 (iGonzalez fc Vanturd[j"998t lOimenezll2000l : 



_ _ — (It . _ _ . 

iHenrv et al.ll2000al: ISantos et all 12004 IChen fe Zhad l200a lEcuvillon et al.ll2006f) . Our theoretical model of 51 Peg 



implies that it may be slightly hotter and more metal-rich than the best-fit parameters obtained by VF05 (see Table [3]). 
The best metallicity estimate for 51 Peg calculated from the theoretical tracks is [Fe/H] = +0.22 . The star is slightly 
older than the Sun, 6.76;i;^ Gyr old. This is in good agreement with FRS99 (7.6 + 4 q Gyr) and most of other isochrone 
analy ses: 4 ± 2 Gyr (jFuhrmann et al.lll998t) . 9.2;^ Gyr ([Nordstrom et alj I2004D and 5.1~ij;g Gyr (jLachaume et all 
1999). From the re-calibrated stellar rotation period determined from the Ca II observation, IHenrv et al.l (12 000a) 
yielded an age of 3 - 7 Gyr. The derived mass, age, radius and the mean rotational period (~ 25 days. IHenrv et al.1 
20004 all indicate that 51 Peg is a star very similar to our Sun, except for its metal-rich atmosphere. 

5.1.3. v Andromedae (HD9826) 

The bright F8V star , i>And is a known trip le-planet system, harboring three Jupiter-size planets at 0.059 AU, 
0.829 AU and 2.53 AU (jButler et al1ll997L Il999). u A nd has a distant sub-stellar companion (M4.5V) at ~ 750 AU 
revealed by co- proper motion (jLowrance et al.ll2002h . which does not affect spectral analysis. Four Lick/Hamilton 



spectra were obtained for v And by VF05. The atmospheric properties derived from each spectrum show a modest 
range of temperature and metallicity, T c ff = 6150 - 6334K and [Fe/H] = 0.08 - 0.192. In general, modeling higher- 
mass stars involve some ambiguity as the sharp hook starts to appear at the end of the main sequence following the 
core hydrogen exhaustion (cf. § 13.21 and Figure fTD|) . The discrepancies in the spectroscopically determined atmospheric 
properties is sensitively reflected in the stellar models, as seen in the derived PDFs. 

Overall, our theoretically derived properties agree with those by FRS99. The derived age-PDF indicates that v And 
is 3.12^Q 2 QGyr old, with two additional probability peaks around it. The majority of the spectral and photometric 
analyses agree with the young age: using theoretical isochrones, 2.9 ± 0.6 Gyr (jLachaume et al.lll999t ). 3.3^'y Gyr 
(iNordstrom et a l. 2004), 3.8 ±1.0 Gyr (jFuhrmann et al.lll998h : using age - activity relation from Ca II flux observation, 



5 Gvr ponaImef t993') ; using the age - [Fe/H] relation, an upper bound of 2.3 Gyr (jSaffe et al.| [2005). Another possible 



model with a relative probability 20 - 40 % is observed in the calculated PDFs, which has a higher mass (~ 1.4M Q ) and 
even younger age (~ 2 Gyr). This is attributed to one of the Hamilton spectra that yields hotter surface temperature 
and higher metallicity than the other three (6334 K and +0.192 ). The SPOCS catalog denotes the mean values of all 
four ob servations, T e g = 6213 K and [Fe/H] = +0.122 , consistent with the other reported v alues, 6107—6212 K and 0.09 
- 0.17 (jFuhrmann et al.l 119981 : 1 Gimenezl 120001 iGonzalez fc Laws! 120001: ISantos et aLll2004l ). With a lower probability, 
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Fig. 9. — Theoretical HR diagram and stellar evolutionary tracks for 51 Pegasi. The tracks correspond to our best estimate metallicity, 
[Fe/H] = +0.22 . The track for O.95M0 is truncated at 14Gyr. The stellar mass and age estimates are 1.05Mq and 6.76 Gyr. The star is 
similar to our Sun except for its metal-rich atmosphere. 
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Fig. 10. — A small subset of theoretical evolutionary tracks used for the modeling of v Andromedae. The temperature and luminosity 
quoted by VF05 are denoted with 1<t error bar and 10ct error box. The tracks represent the evolutionary sequences of stars with metallicity 
[Fe/H] = +0.20, corresponding to our best-estimate metallicity. The star is most likely a 1.31A/q main-sequence star that is 3.12 Gyr old. 
Another model with M < 1.25Mq at a main-sequence turn-off age (~ 5 Gyr) is also possible, although with a much smaller probability. 

v And could also be a near turn-off star with a lower mass (~ 1.25M©) and an older age (~ 4.5 Gyr), not because of the 
observational ambiguity but because of the sharp rise in the stellar luminosity after core hydrogen exhaustion. This 
probability is relatively smaller due to the rapid evolution of the star away from the main sequence toward the lower 
temperature. As mentioned above, all other isochrone analyses with appropriate statistical treatment (e.g., Bayesian 
analysis) favor the main-sequence model younger than 3.5 Gyr, rather than turn-off or post-main-sequence models. 

5.1.4. 55 Cancn (HD75732) 

55 Cnc is the richest p l anet-host star known t o date, harboring four Neptune- to Jupiter-size planets (|Butler et al.l 
fl997t lMarcv"et~aT][200a IMc Arthur et~all[200l . This is a system of great interests for many aspects of planetary 
dynamics. The innermost planet 55 Cnc e is a hot sub-Neptune mass planet, while the outermost planet 55 Cnc d is 
one of the few planets known with an orbital se mimajor axis comparable to that of Jupiter (a = 5.257 AU). The middle 
two planets are in 3:1 mean motion resonance (jJi et al.ll200llMarzari et al.ll2005h . posing an interesting question for 
dynamical stabilit y. 55 Cnc is also a visual binary system, with a stellar companion at ~ 1100 AU away from the 
planet-host star : I lollleii fc Jascheklll982t iMugrauer et al.ll2006f). In addition, the star was onc e claimed to have a 
Vega-like dust disk based on infrared observations |Dominik et al.l Il998t iTrilling fc Brownl |1998[ ) . Observations with 
different wavelengths have set an upper limit on the disk mass much lower than the previous estimate from the 
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Fig. 11. — Observed temperature and luminosity of 55 Cancri and the evolutionary tracks of stars with the best-estimate metallicity 
[Fe/H] = +0.38. All the tracks are terminated at 14 Gyr. Effective temperatures determined from other spectroscopic observations typically 
lie within a range, log (T e f[/K) = 3.70 — 3.73. The star is most likely a main-sequence star with a mass 0.967q'q5Mq 



infrare d observation, and now it is generally agreed th at the infrared excess most likely came from a background 
object (jJavawardhana et alj|2002t iSchneider et al.ll2001[ ). 

The strong photospheric absorption lines of 55 Cnc indicate an anomalous metal abundance, classifying the star as 
a so-called "super-metal-rich" star. The peculi ar spectrum of 55 Cn c results i n controversia l atmospheric properties. 
The star is normally classified as a G8V star (jCowlev et al.l fl967'). whereas iTavlorl (|1970t ) identified it as a super- 
metal-rich K dwarf. Rep orted surface temperature and metallicity are not yet well constrained: T e g = 5100 — 53401C 
and [Fe/H] = .20- 0.45 (lArribas fc Martinez RogeHll989tlBaliunas et alJ ll997t iFuhrmann etlHll998t lGonzaledll998l 
lGimenez|[2000l: iReidl l2002t iSantos et al.ll2004t lEcuvillon et all 120061 ). Many models show that 55 Cnc is a sub-solar 
mass main-se quence star, whereas some other observations claim i t to be a subgiant, based on the atmospheric CN 
enhancement dGreenstein fc Oinasl 119681 : lTavlor|[l970t lOinasI [1971) and relatively low s urface gravity, logq = 4.10 
()Gonzalezlll997| ). However, using the greatly increased number of Fe I and Fe II lines, iGonzalez fe Vanturel (|1998f ) 
revised the surface gravity to loggr = 4.40, which is too large for a normal subgiant star. For our models, we have 
used the nine Lick/Hamilton spectra obtained by VF05. The mean temperature and metallicity yielded by SME are 
T off = 5253 ± 44 K and [Fe/H] = +0.31 ± 0.029 . 

The reported physical parameters of 55 Cnc also show discrepancies. Most of the analyses claim a model that is 
either very young and slightly more massive than the Sun (~ 1.05M Q ) or very old and slightly less massive than 
the Sun (~ O.9OM ). Some of the spectroscopic analyses yield an age close to the Hubble time (|Perrin et al.l 119771 : 
ICavrel de Strobell Il987t IGonzalez! [1998) . FRS99 also predicts an old age 8_|_7 Gyr and a sub-solar mass O.95M0. 
Naively speaking, an extremely old age of 55 Cnc seems unreasonable considering its anomal ously metal-abundan t 
atmosphere. Chromospheric Ca II activity of 55 Cnc suggests inste ad a young age of 5 Gy r (|Baliunas et al.|[l997f) . 
Using the high metallicity +0.40 and keeping the solar He-abundance. IFuhrmann et al.l (|1998l ) derived an age < 5 Gyr. 
It has even been su ggested that 55 Cnc might be a member of the Hyades Supercluster with an age upper limit of 2 Gyr 
(|Eggenlll985lH992l ). However, various observational clues such as the measured stellar rotational period of 41.7 days, 
possibilit y of an extreme H-deficien t atmosphere and the lack of detectable lithium, argue against extreme youth of 
the star (|Gonzalez fe Vanturdll998l ). The substantial size of the convective zone also supports a more evolved stellar 
model. 

Among the five sample stars we discuss here, 55 Cnc has the least-constrained mass and age. The observations by 
VF05 suggest it is a main-sequence star with M — O.96^g'Q5M and [Fe/H] = +0.38 (see FigurefTTj). Due mainly to the 
slow main-sequence evolution of metal- rich stars with M < l.OM , the derived age of the star is poorly constrained. 
Although the age-PDF is indicative of an older age, nearly any age is po ssible within the r ange of 0-14 Gyr (Figure[7]). 
A choice of higher temperature and metallicity (e.g., 5336 K and +0.40 . IFuhrmann et all ([T998)) would favor a younger 
age and a larger mass. The situation can be understood as a result of the assumption for He-abundance. The scaling 
of He with respect to the other metals (SY/6Z) is not well known, especially for these super-metal-ric h stars. Using the 
stellar evolutionary tracks with the high metallicity but keeping the solar He-abundance (SY/SZ = 0) . IFuhrmann et all 
(|1998h sho wed that 55 Cnc is located b elow the ZAMS line in the HR diagram. If the H-depletion in the atmosphere is 
confirmed (|Gonzalez fc Vanturel 1998), He-fraction needs to be scaled up as SY/5Z = 2.5 ± 1.0, as adopted by FRS99 
and us. Using a uniform prior distribution of 5Y/6Z, our best fit for 55 Cnc occurs for our tracks with 5Y/SZ = 1.5. 

The model posterior temperature and metallici ty, 5327 K and + .38 a re considerably higher than the values from the 
SPOCS catalog but similar to the observation bv lFuhrmann et all |l998). The derived surface gravity log g — 4.48^Q'gg 
is quite large and in agreement with VF05 and FRS99. 
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Fig. 12. — Theoretical HR diagram and stellar evolutionary tracks for pCoronae Borealis. Tracks with our best-fit metallicity, [Fe/H] = 
+0.17, are shown. The tracks with M = 0.90Mq and 0.95-Mq are truncated at the maximum age of the grids, 14Gyr. The derived model 
of pCrB yields M = 0.96 ± 0.02A/© and r = 11.04~g;|| Gyr. 

5.1.5. p Coronae Borealis (HD143761) 

The AFOE (Advanced Fiber Optic Echelle spec trograph) team discovered a Jupiter mass planet in 39.8 days orbit 
around the GO - 2 V star pCrB (jNoves et al.lll997aH bl). It is the only star with sub-solar metallicity among the five stars 
discussed in this section. The mean metallicities and effective temperature of the star yielded by one Keck/HIRES 
spectrum and three Lick/Hamilton spectra by VF05 are T c s — 5823 K and [Fc/H] = —0.14. Other spectroscopic 
analyses all identify the metal-poor atmosphere of pCrB, ranging from [Fe/H] = —0-19 to —0-32 dGratton . et all 
1996t iKunzli et al l|1997t IFuhrmann et al.l ll998l; lGonzalez|[T998t lGirnenezll2000j: iHenrv et al]|2000at iTakeda et al.ll2001l ; 
Santos et al.ll2004t lEcuvillon et al.ll2006l ). 

FRS99 yielded a model with an extremely old age, 14 ± 2 Gyr, and M = 0.89 ± 0.05Mq. Our analysis predicts a 
model that is slightly younger and more massive, 11.04^^ Gyr and 0.96 ± 0.02M Q . This mass is consistent with 
the isochrone mass of 0.95M Q by ISantos et al.l ([2004) , derived from the SARG spectrum. The ne ar turn-off main- 
sequence age has been confirmed by other isochrone / evolutionary track analyses; 10.2 ± 1.7 Gyr (iFuhrman n et al.l 
11951 ): 12.1"!° Gyr (jNordstrom et al.ll200l ; 12.1 ± 0.9 Gyr (|Ng fc Bertellil fl998l) . The paucity of heavy elements 
detected i n the atmosphere is consistent with some of these very old ages, although metallicity is typically a poor age 
indicator (Saffe et al. 2005). Other observational evidences support that pCrB is an evolved, near sol ar-mass star. 
Relativ e ly low chromospheri c activity and slow r otational period (~ 20 days) have been confirmed by iNoves et alJ 
(|1997af ): [Henry et"aT1 (l2000al). INoves et a l. (1997a) suggests that the high proper motion of pCrB out of the Galactic 
plane at 28km/s (jCavrel de Strobe] 1 1996| ) may indicate that the star was a member of the old disk population. Our 
model with the slightly developed convective zone also supports that the star is at least near the end of the main- 
sequence stage. 

The theoretically derived atmospheric parameters are consistent with the spectroscopic values by VF05, except that 
the model posterior surface gravity is significantly lower, lo g. 17 = 4.18. Th i s is co nsistent with the is ochrone value of 
4.14 b y VF05, and other spectroscopic observations, 4.11 bv lGratton et al.l (|1996l ) and 4.05 - 4.19 bv IFuhrmann et all 
lfl998h . 

5.2. Sample Models for Other Known Planetary Systems 

In addition to the stars previously modeled by FRS99, here we present the models for 5 planet-host stars with 
particularly interesting stellar or planetary properties that are worth detailed discussions. Table[5]lists the spectroscopic 
parameters and the derived physical properties of these stars. 

5.2.1. HD209458 

HD209458 is the first extrasolar planetary system for which transit events were observed (jCharbonneau et alJl2000t 
IHenrv et~a l. 2000b). Accurate stellar modeling for stars with transiting planets is crucial for determining the radius and 
hence the interior structure of the planet. The tightly constrained orbital inclination angles achievable for transiting 
planets can eliminate the factor of sini from the observed radial velocity and yield an exact planetary mass (M p i). 
From the photometric curve and the theoretically determined stellar radius, the radius and therefore the density (and 
c omposition) o f the p lanet can be found. 

iMazeh et al.l (|2000t ) have created several stellar models using isochrones derived from the Geneva, Padova, Claret 
and Yale code. By locating the observed My and T e ff and interpolating between the isochrones, they determined 
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TABLE 5 

Sample Models for Stars with Known Planetary Companions 





Observed Data 




Posterior Model Parameters 




Star 


T cff [K] [Fe/H] logg [cgs] 


M„[M ] 


Age [Gyr] i?,[R ] M ce [M ] flce[R ] a 


log g [cgs] 



HD177830 


4949 


+0.33 


4.03 


HD209458 


6099 


+0.02 


4.38 


HD27442 


4846 


+0.29 


3.78 


HD38529 


5697 


+0.27 


4.05 


HD69830 


5361 


-0.08 


4.46 



> 1.91 3.24 



1.13 
1.59 



-0.02 
+0.03 
-0.14 
=+0.09 

1.48+0.05 
0.85+0.01 



+().-)(> 
2 44" 1S4 

2 - 84 +0.60 



3.28 



0.24 
+0.36 
> 12.04 



2.95~ 
1.14 
3.43 

2 - 50 +o.os 
0.90+0.02 



0.945] 



+o.O() 

0.03 
+0.11 
0.06 



-0.021 
J +0.114 

0.007+0.002 
1 014-°' 063 

l.U14 + () ()48 

064~ (M,(,U 
U ' UD ^+0.019 

0.038 +0 002 



> 1.29 

24" 01 
u -^*+0.02 

> 1.87 
7i-° 01 

u - ' 1 +0.03 

0.29+0.01 



3.91+0.01 
4.39+0.04 
3.56] 



-0.04 
-+0.05 

3.94+0.02 

-0.01 



4.47 



+ 0.02 



Note. — Theoretical models for stars with known planetary companions discussed in § 15.21 Adopted uncertainties for the observed 
parameters are, 44 K for T e ff, 0.029 dex for [Fe/H] and 0.060 cm s -2 for logg. 
a Depth of the convective zone, measured from the outermost stellar surface 

M = 1.1 ± CUM© and R = 1.2 ± 0.1i? o . We have derived a similar model with better constraints, M - 1 t3-° 02 



and R = 1.1 4: + ^^Rq. Using th e stellar parameters derived by iMazeh et al.l (|2000f ) and the observed inclination i = 
87^1±0 o .2. ICharbonneau et all (pOOOl) determined the planetary parameters, M n \ = 0.63M , T and i? pl = 1.27±0.02i?j. 
Using the theoretical isochrones bv lBertelli et alj (ll994l).lAllende Prieto fe; Lambert (|1999t l determined a stellar radius 
which is closer to our result, R = 1.15 ± 0.08-Rq. iHenrv et al.1 1 2000bD have adopted this model and derived a slightly 
larger planetary radius R p \ = 1.42 ± 0.10i?j than the value bv lCharbonneau et alj (|2000f ) . 



5.2.2. HD 69830 

The Doppler detection of three planets around the star HD69830 was recently reported by lLovis et alj ((2006) . 



It 



is the first triple-planet system consisting only of Neptune-mass planets: M pl sini = 10.2Af ffi (planet b), 11.7 
(planet c) and 18. 1M® (planet d). Prior to the planet detectio n, the system attracted a lot of interest owing to the 
large infrared excess observed by the Spitzer Space Telescope (|Beichman et alj|2005al) . indicating the presence of a 
massive asteroid belt within 1 AU. 

Using the high resolution spectra obtained with HARPS at La Silla Observatory. lLovis et al.l (|2006( l determined an 
effective temperature T cf f = 5385 ± 20 K and metallicity [Fe/H] = -0.05±0.02. The quoted values of T eff = 5361±44K 
and [Fc/H] = —0.08 ± 0.03 by VF05 are in good a gree ment with thes e value s. The isochrone analysis using the 
theoretical evolution models by iSchaller et al.l (|1992h and iGirardi et al.l (|2000h yielded a stellar model of HD69830 
with a mass 0.86 ± 0.03M Q a nd an age ~ 4 - 10 G yr. Our calculation shows a similar model with a mass 0.85 ± 0.01Af Q 
and an older age > 12 Gyr. lLovis et alj (|2006f ) also performed numerical TV-body simulations to test the dynamical 
stability of the system assuming coplanarity of the orbits. They tested two inclination angles i = 1° and 90°, and 
in both cases the system remained stable for at least 1 Gyr. Note that long-term stability lasting for as long as our 
derived lower-limit age 12 Gyr might favor smaller planetary masses, corresponding to a near edge-on view of the 
system (i ~ 90°). A more extensive stability analysis as well as a better-constrained stellar age are needed to provide 
tighter constraints on the orbital properties of the planets in this system. 

5.2.3. HD27442, HD38529 and HD177830 

The SPOCS catalog contains 86 well-observed subgiants, and about 10 of these subgiants have detected planets. 
Interestingly three of them are observed to be super-metal-rich subgiants. We find it worthwhile to present our 
theoretical models for these subgiant stars because their physical properties are usually not well-determined. It is 
particularly difficult to derive an accurate stellar mass for subgiants, since in the HR diagram theoretical subgiant 
tracks with different masses are much more closely separated than main-sequence tracks (see Figure [1}. On the other 
hand, derived ages of subgiants are relatively well-constrained, because of the rapid cooling of the stellar atmosphere 
toward the r ed-giant phase. 

HD27442 (|Butler et al.ll200l is one of the most evolved stars in our stellar sample, possibly already at the beginning 
of the red-giant phase. It is a K 2 IV a star with a planet HD27442b with minimum mass M p \ sin i — 1.35Mj and period 
P = 423.8 days. The observed parameters of t he stars from the SPO CS catalog are T cS = 4846 K, [Fe/H] = +0.29 and 
log g = 3.78. From spectroscopic observations, Randich et al. (1999) derived the atmospheric parameters T e ff = 4749 K, 
[Fe/H] = +0. 2 2 and \ogg = 3.3. They further combined their observations with theoretical isochrones computed by 
iBertelli et al.l (|1994l ). and determined the mass 1.2 ± O.IMq and age r = 10 Gyr by interpolating between two sets of 
isochrones, [Fe/H] = 0.00 and +0.40. We have derived a model that is younger and more massive than their model, 



M 



1.48;°;^M Q and r 



2.84 



-0.36 
+0.60 



Gyr. 



HD38529 is a G4IV sub giant star harboring two giant planets. The inner planet HD38529b (M pl = 0.78Mj, P = 
14.3 days) was discovered bv lFischer et al.l (|2001l ). A long-period residual trend reported at the time was l ater confirmed 
to be another planet HD38529c, with a mass M p \ = 12.7Mj an d an orbital period P = 2 1 74 da ys (jFischer et al.l 
120031 ). The planetary masses are derived using the stellar model bv lAllende Prieto fc LambertJ (| 1999( 1: M* = 1.39M , 
i?» = 2.82i? Q and logg = 4.13. Interestingly, the large mass of HD38529c is close to the theoretical deuterium burning 
limit and thus suggests that it may be a sub-stellar companion rather than planetary. We have derived the stellar 
mass 1.48 ± 0.05M Q , indicating a companion mass even larger than 12.7Mj. 
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Fig. 13. — Age vs mass for 669 stars with well-defined ages. The 74 planet-host stars are marked with red open circles. 



The KOIV star HD 177830 is a highly evolved subgiant with very high metallicity, [Fe/H] = +0.33. A J upite r-size 
planet HD177830b (M pl sini = 1.22Mj) with an or bital period of 391.6 days was discovered bv lVogt et all (|200Gf ). By 
interpolating the theoretical evolutionary tracks bv lFuhrmann et all (|1997lfl998l ) to the observed My and B — V, they 
estimated the stellar mass of 1.15 ± 0.2Af Q . However, a mass of a subgiant star estimated by interpolating between 
evolutionary tracks can often be inaccurate. Moreover, modeling a star with such high metallicity requires a grid of 
theoretical evolutionary tracks with a high metallicity resolution and a range of helium fractions. We have derived 
a mass M > 1.91M© and a radius R = 2.62;^i? . There is a small local maximum at 1.55Af© in the mass-PDF 
with a ~ 10% probability relative to the global maximum at 2.0M Q . The age of HD177830 is well constrained: 
r = 3.24^Q 5g Gyr. HD177830 is a highly evolved super-metal-rich subgiant, located very closely to HD27442 in the 
HR diagram. However, HD177830 is even more metal-rich ([Fc/H] = +0.33) than HD27442, thus its mass needs to 
be at least comparable to that of HD27442 or even larger in order to reach the same evolutionary stage as HD27442 
within a similar but slightly older age of ~ 3 Gyr. 

5.3. Parameter Correlations 

In this section we explore possible correlations between the derived stellar parameters. It should be remembered 
that there are many spurious trends that are artificially introduced during the stellar modeling process or because 
of observational selection effects. Those artifacts need to be carefully removed to analyze any meaningful statistical 
correlations between the derived parameters. 

5.3.1. Age - Mass Relations 

Figure [TBI shows age - mass relations for the SPOCS sample. Only the 669 stars with well-defined mass and age are 
plotted, to remove the artificial accumulations at extreme ages (cf. § !3.4.ip . 

Two major features are discernible in the figure: a large envelope of main-sequence dwarf stars extending from the 
upper- left (young, high-mass stars) to the lower-right corner of the figure (old, low-mass stars), and ~ 14 subgiant 
stars, in the region M > 1.4M Q and r > 4 Gyr. The shape of the large envelope is mostly determined by observational 
selection effects and our stellar modeling method: (i) The upper age limit for each stellar mass roughly corresponds 
to the main-sequence lifetime of the given stellar mass (e.g., ~ 10 Gyr for 1M@). As the stellar mass increases, the 
age of main-sequence stars is typically better constrained. Note that the wide range of metallicities creates a scatter 
in the main-sequence lifetime, since metal-poor stars typically evolve more rapidly (e.g., a IMq star with only 10% 
of solar metallicity leaves the main sequence after ~ 7 Gyr). (ii) Each age has an approximate lower boundary with a 
critical mass below which well-defined age cannot be determined. In general, for low-mass, unevolved stars near the 
ZAMS, only an upper bound of age can be derived, since these stars slowly evolve upward in the HR diagram and 
do not leave the main-sequence track within 14 Gyr. A few exceptions such as HD144253 and HD65583 in the region 
M < O.85M and r < 4 Gyr have particularly low metallicities, [Fe/H] = -0.21 and -0.48, respectively. Although the 
derived la age uncertainties for these two stars are still large (8.7 Gyr and 8.2 Gyr, respectively), these metal-poor 
stars experience a more rapid rise in luminosity during the main-sequence phase which helps constrain the ages better 
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Fig. 14. — Derived well-defined age vs metallicity. Here, red open circles represent 74 planet-host stars with well-defined ages. An 
integrated age-PDF is calculated from all the age posterior PDFs for each metallicity bin, in order to estimate the age dispersion for given 
metallicity. The mode value and 68 % credible interval of the stellar age in each bin are shown for the entire sample (solid square) and for 
the planet-host star sample (open square). 



than for slowly evolving metal-rich stars, (iii) The narrow void on the left of the envelope is a lack of well-defined ages 
less than ~ 1 Gyr, mainly caused by the blue color cutoff in the SPOCS catalog (B — V > 0.5) and exclusion of stars 
that are chromospherically very active. 

Similar trends are seen in the sample of planet-host stars. Planet-host stars are dominant in the mass range 
M = 0.95 — 1.4M . Doppler radial-velocity observations become more challenging for stars with masses above this 
range because of the increasing atmospheric jitter. Planetary companions are not common around low-mass later-type 
stars (M < O.9M ), although current spectroscopic surveys can achieve precise radial velocities for such low-mass 
stars. 

5.3.2. Age - Metallicity Relations 

It has been confirmed by various spect r oscopic observations that planet-host s t ars are on average more metal-rich 
than single field dwarfs (iGonzaled Il997t iFuhrmann et al.l Il997t lGonzale3 119981 : ISantos et al.l 120001 : I Gonzalez et al.l 
1200 It ISantos et al.l 120011 l2003f ). One explanation posits that planets are more efficiently for med around stars with 
higher metallicity because of the higher fraction of solids available in the circumstellar disk (jlda fc Linll2004bh . An 
alternative theory suggests that the observed high metallicities of planet-host stars are caused by late-stage accretion 
of gas-depleted material. In this "pollution" or "planet-accretion" hypothesis, solid bodies (e.g., planetesimals or giant 
planet cores) migrate into the stellar atmosphere and enhance the stellar surface metallicity. 

Theoretical calculations show that the degree of observable metal-enhancement is largely dependent on the size of the 
stellar convective zone (and therefore the effective temperature), but independent of the stellar age. ICodv fc Sasselovl 
(2005) tested the effect of planet accretion on the subsequent stellar evolution by calculating various stellar evolution 
models with polluted stellar atmospheres. They added metals to stellar convection zones at various arbitrary times 
up to ~ 6 Gyr and showed that a polluted star will reach the same equilibrium state, regardless of whether the 
metal-rich material is accreted immediately or at a later time. Their calculations suggest that polluted and thus 
metal-enhanced stellar atmospheres cannot be distinguished from intrinsically metal-rich stellar atmospheres in age - 
metallicity relations. 

Figure Q3] shows an age - metallicity scatter plot for 669 stars with well-defined ages, including 74 planet-host stars. 
It should be remembered that many planet-search programs preferentially select metal-rich stars to optimize planet 
detections, thus our sample is not expected to represent the true age - metallicity relation for solar-neighborhood 
stars. Another caveat is that age is the most poorly-determined posterior parameter, thus a simple scatter plot for 
the mode ages and observed metallicities is not an accurate indicator of any possible age - metallicity relation. To 
see qualitatively the distributions of the derived ages and the observed metallicity, we have computed a sum of the 
normalized age-PDFs for each given metallicity. In Figure HH the mode and associated la uncertainties for the 
integrated age-PDF is presented for each metallicity bin. Within the metallicity range [Fe/H] =-0.3 - 0.4, the age is 
nearly uniformly scattered around the solar age (~ 4.5 Gyr) in each metallicity bin. Note that the apparent mode-age 
- metallicity relation is not a real trend but is an artifact of the age - mass relation (Figure [T3|) and the mass - 
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Fig. 15. — The mass and width of the convective zone vs metallicity. Here R CE is defined as the distance between the stellar surface and 
the bottom of the convective zone. Planet-host stars are marked with red circles. The apparent deficit of stars with M ce > O.IMq comes 
from the relative shortage of K-type dwarfs in the SPOCS catalog and the steep decline of the convective envelope mass for stars earlier 
than F-type. A possible weak correlation is present between M ce and maximum [Fe/H] in the range 



metallicity relation (Figure IT5|) . 

While we do not find a significant age - metallicity relation in the planet-host star sample, this may be largely 
due to the limitations of our analysis technique. Several other studies of the age - metallicity relation for planet-host 
stars also showed li ttle or no variation in metallicity as a function of age ijSaffe et al.l 12005; Bci chman et al.l [2005b; 
iKarataset aLl l2005V The mean metallicity of the planet-host stars in the sampl e is hig her by 0.12 than that of all the 
stars with well-defined ages, consistent with the analysis by iFischer &; Valentil (|2005l ) Interestingly, there are three 
planet-host stars with old ages (r > HGyr) and extremely high metallicity: HD30177 ([Fe/H]=+0.48), HD45350 
([Fe/H]=+0.30) and HD73526 ([Fe/H]=+0.27). These three stars are all super- metal- rich late G-type stars with near 
solar mass, very similar to 55 Cnc. The fact that most of the planet-host stars with very old ages (r > 10 Gyr) have 
super-solar metallicity largely co mes from observational selection effects: (i) the frequency of planetary companions 
increases with stellar metallicity (jFischer fc V alcnti 2005|) , and (ii) stars with sub-solar metallicity evolve more rapidly. 
Metal-poor F-,G- and K-stars most likely already left the main sequence before 10 Gyr and thus usually are excluded 
from planet-search programs. 

5.3.3. Metallicity - Convective Zone Relations 

The derived mass and depth of the convective zone are plotted against the observed metallicities in Figure [151 The 
noticeable scarcity of stars with M ce > 0.05M Q most likely comes from the relative s h ortag e of K-type stars in the 
sample, corresponding to T e g < 5500-fT (see Figure 12 of VF05). iPinsonneault et all ((2001) computed the mass of 
convective zone for stars with masses in the range 0.6 - 1.3M Q and showed that M ce is a sensitive function of the 
effective temperature. For example, they showed that the convective envelopes of F-stars are more than 10 times less 
massive than that of K-stars (also see Figure [16]). The sharp decline of M ce for stars earlier than K-type, combined 
with the relative shortage of K-stars in the SPOCS catalog emphasizes the decrease in the population of stars with 
M cc < 0.05M© in Figure [151 

The distribution of M cc can be also used to test the metal-enrichment mechanism by planet accretion. If planet 
accretion is a common phenomenon among planetary systems, we would expect an increase of maximum metallicities 
toward lower values of M ce because mixing of accreted material is less effective in a thinner convective envelope. Thus, 
the enhanced metallicity is better preserved for stars with less massive convective envelope. In particular, if planet 
accretion is responsible for the high metallicities of planet-host stars, then K-dwarfs and subgiant stars should show 
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Fig. 16. — Stellar mass vs metallicity (left) and stellar mass vs convective envelope mass (right) for the entire sample. The planet-host 
stars are marked with red circles. Note that M ce starts to decrease faster for stars with masses above I.OMq. 



systematically lower metallicities, since the large convective envelopes of these stars completely dilute the accreted 
solids. 

Figure [15] shows a weak positive correlation between the maximum [Fe/H] and Af ce , which is opposite of what 
is expected. The most metal-rich planet-host stars in the sample are not F-type stars, but K- or G-dwarfs with 
larger convective zones: HD145675 (KOV, [Fe/H] = +0.41), HD2039 (G2/3IV - V, [Fe/H] = +0.39) and HD30177 
(G8V, [Fc/H] = +0.34). Note that F-type stars in the sample (M — 1.2 — 1.5Af ) are preferentially more metal- 
rich, because metal-poor F-stars have shorter main-sequence lifetimes (see Figure fl6|). Despite this selection bias, 
the maximum metallicity and mean metallicities of the F-type stars in the sample are lower than those of the G- 
and K-type stars. Furthermore, all three planet-host subgiants with well-defined ages are in fact super-metal-rich: 
HD27442 ([Fe/H] = +0.39), HD38529 ([Fe/H] = +0.31) and HD177830 ([Fe/H] = +0.36). This cannot be attributed 
to planet-pollution mechanism because subgiant stars develop large convective envelope in which accreted materials 
are well mixed with the deeper stellar interior. Although a larger sample of K-dwarfs and subgiants would be more 
desirable, these results argue against significantly greater (> 0.1 dex) pollution of the convective envelope of stars with 
planets than stars without planets due to accretion of planets or planetesimals. 

5.3.4. Stellar Mass - Planetary Mass Relations 

In spite of the chaotic nature of the formation and subsequent dynamical evolutions of planetary systems, it has been 
suggested that the distributions o f planetary param et ers ma y be mostly controlled by several key stellar parameters. 
Extensive work has been done bv llda fc~L in (2004a. b, l2005f) in search for a deterministic theory of planet formation. 
Planet-formation models can also predict subsequent orbital evolutions of planets such as migration or orbital decay 
by tidal dissipation, based on the initial configurations of the protoplanets. Later evolution of planetary orbits may 
also be characterized by the stellar properties. Thus, the derived distributions of the planetary and stellar properties 
might indicate correlations that constrain certain formation and dynamical theories. 

Figure [171 shows the relation between the observed planet mass (M p isini) and the derived stellar mass (M„). 
The current planet formation theories based on planetesimal coagulation strongly depend on the properties of the 
circumstellar disk. The growth of planetesimals is highly dependent on the surface density of the disk, thus stars 
with initially more massive disks are expected to produce larger planet masses. However, the distribution of the disk 
masses inferred from infrared observations of T-tauri stars is not yet well-constrained. The observations of dust in 
protoplanetary disks indica te that the total disk mass of a solar-type star can typically range from 10~ 4 to 1O _1 M 
(|Beckwith fc SargentJll996T) . Thus, accurate mass of protoplanetary disks cannot be determined as a function of stellar 
mass for the stars of our interest (M = 0. 5 - 2.0M ff ) . 

Using extensive numerical simulations, llda fc LinT (|2005[ ) derived a positive correlation between the characteristic 
planetary mass and the stellar mass within the range 0.2 - 1.5Af Q . In Figure [T71 the maximum planetary masses 
do increase in the same range of derived stellar masses. However, the sample of only five systems is far too small 
to assert any correlation. Also note that there is an intrinsically different population of companions above ~ 13A/j, 
corresponding to brown dwarfs. Spectroscopic observations show that there is a distinct absence of secondaries with 
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Fig. 17. — Stellar mass vs planetary mass. For multiple-planet systems, the combined masses of all the planets are marked with circles. 
The size of the circles increases in the order of double-, triple- and quadruple-planet systems. The shaded region corresponds to the "brown 
dwarf desert". The two systems above M p isini = 19Mj, HD202206 and HD168443, are probably triple systems with a planet and a brown 
dwarf companion. The companion of the star HD162020 with a minimum mass M p isini = 14AM] is also likely a brown dwarf. The four 
stars with masses > 1.6Mq are evolved subgiants. 



masses M 2 0.01 — 0.08M (n around solar-type primaries, so-called "brown dwarf desert" (Halbwa chs et al] 120001 : 
iGrether fc Lineweaverl [20061 ) . This suggests a completely different formation channel for planets than for brown 
dwarfs. Interestingly, the three systems with total companion mass beyond 13Mj in Figure [l"7l HD162020, HD168443 
and HD202206, may be the rare candidates in which a K- or G-dwarf primary contains a brown dwarf companion. The 
measured compan ion masses are, M pl sini = 14.4Mj (HD162020b), 17.5Mj (HD202206b) and 16.9Mj (HD168443c) 
(|Udrv et al.ll2002t ). HD202206 and HD168443 are triple systems, associated with another massive planet, HD202206c 
(2.44Mj) and HD168443b (7.7Mj). It has not yet been confirmed whether these massive companions are sub-stellar 
or planetary. HD162020b is particularly interesting, as it has a tight orbit similar to hot Jupiters (P ~ 8 days). If it 
is confirmed to be a planet, that would significantly change the upper bound of the M* - M p isini distribution. 

Lastly, Figure [17] shows a lack of massive planets around massive stars (M* > I.8M0). Doppler planet detection is 
intrinsically more challenging for stars earlier than F-type, and thus the sample in this region of parameter space is 
fairly incomplete. Nevertheless, radial-velocit y searches should fayo r detections of heavi er planets. Currently ongoing 
plane t search programs targeting A - F stars (Galland et al. 2005b a|) and G - K giants (|Sato et al.ll2005t lHatzes et al.l 
|2005( ) will fill more samples in this region and may reveal interesting trends. 

5.3.5. Stellar Properties and Planetary Orbits 

Planetary orbits can still drastically evolve after the early phase of planet formation. As a result, the distributions 
of orbital parameters reflect a mixture of various long-term dynamical histories as well as the initial configuration of 
planets, which may be partly obscured. Thus, it is quite challenging to interpret any possible relations between orbital 
a nd stellar properti es. 

lUdry et'a l. (2003) pointed out a lack of massive planets (M p isini > 2Mj) with short periods (P < 100 days) in the 
observed sample of extrasolar planets. Interestingly, a similar void may be observed in our stellar mass - orbital period 
relation. Figure [TBI shows a deficit of short-period planets around massive stars. One possible explanation for the lack 
of high-mass close-in p lanets is that Type II migration is less effective for more massive planets because of a larger 
gap-opening timescale ([Trilling et al.ll2002| ). An alternative mechanism to eliminate short-period m assive planets is th e 
orbital decay due to tidal dissipation, although this is only effective for orbital periods < 10 days ([Rasio et al]|1996l) . 
Note that the five most massive stars above the void in Figure [18] are all subgiants. As a star evolves toward the 
red-gia nt phase, it develop s a deep, massive convective envelope which more effectively dissipates the planet's orbital 
energy. iRasio et al] (|1996f) demonstrated that if the star 51 Peg expands to about twice the Sun's radius, the orbital 
decay timescale drastically decreases by a factor of ~ 10 4 , and inevitably the planet will be swallowed in the stellar 
envelope. Since the orbital decay time is also sensitively dependent on the initial orbital semimajor axis, planets at 
initially smaller distances have smaller chances of survival. Since tidal dissipation is effective only within the range of 
a few stellar radii, it cannot explain the lack of planets with intermediate orbital periods (P = 10 - 100 days) around 
massive stars. Although the larger probability of fast Type II migration during the formation and later tidal orbital 
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Fig. 18. — Stellar mass vs orbital period of planets. Multiple-planet systems are denoted by solid circles. Planets orbiti ng around a 
compo nent of a stellar binary are denoted by open stars. The period of 100 days corresponds to the boundary observed by lUdrv et al.l 
(2003) below which no massive planet (M p \ > 2Mj) is detected. 




Fig. 19. — Derived stellar ages vs orbital eccentricities of planets. Planets in multiple-planet systems are marked in solid circles. The 
open stars represent planets orbiting around a component of a stellar binary. Age uncertainties (68 % credible interval) are shown for 
planets with eccentricities larger than 0.6. There is a lack of planets with very high eccentricity (e > 0.7) and young age (t < 4Gyr). 



decay may partly explain the lack of short-period planets around massive stars as well as the lack of short-period 
massive planets, more detections of planets around young massive stars and evolved subgiants are needed to further 
explore these ideas. 

Figure [TO1 shows the orbital eccentricities of 120 extrasolar planets (99 systems) in the sample and the derived ages of 
the systems. One of the striking orbital characteristics of extrasolar planets is their high orbital eccentricities. As of May 
2006, the mean eccentricity for 175 known extrasolar planets is 0.25, higher than that of Mercury (0.21) or Pluto (0.24) 
in our Solar System. Planets formed in the standard scenario are expected to have a nearly circular orbits, thus the 
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high eccentricities of extrasolar planets require later eccentricity perturbation mechanisms (e.g.. lTremaine fc Zakamskal 
2004!, and references therein). 

Although there is no apparent age - eccentricity correlation observed in Figure 1191 there are a few signatures in 
the age - e distribution that may indicate eccentricity evolution of the planets, (i) Most of the very eccentric planets 
(e > 0.6) are in single-planet systems. Two such planets are orbiting a component of a known stellar binary (HD80606 
and 16 Cyg). (ii) Most of the very high eccentricities (e > 0.6) are observed in old systems (r > 5Gyr). (iii) Multiple- 
planet systems have a wide range of eccentricities with an upper bound around e = 0.6. Systems w ith multiple giant 
planets of comparable masses c an easily achieve high eccentricities through dynamical instabilities (|Ford et al.ll2005f ) 
or crossing orbital resonances (|Lee fc Pealdl200l IKlev et al.ll200l . However, very high eccentricities can also lead 
to orbital crossings, and the planets ending up colliding or being ejected from the system. While a larger sample 
is clearly needed, we find this particularly interesting since the planet-planet scattering model (with two planets 
on init ially circular orb its) predicts a rapid decline in the frequency of eccentricities above ~ 0.6 and none above 
~ 0.8 (jFord et al.l 120031 ). Secular perturbations from distant companions can also evolve planetary orbits into very 
high eccentricities. When a planet orbits around a component of a wide binary system with a sufficiently large 
relative orbital inclination (i > 40°), the planet's orbit undergoes a long-term eccentric ity oscillation with a maximum 
eccentricity up to almost unity {Kozai 19621: llnnanen et"aT]|1997tlTakeda fc R asio 2005). Theoretical models show that 
this secular perturb ation (the "Kozai mecha nism") can likely ex plain the high eccentricities of HD80606b (e = 0.93, 
IWu fc Murray! [2003T ) and 16 Cygb (e = 0.67. lHolman et aT]|1997f) . Although the period of the eccentricity oscillation 
is dependent on the semimajor axis and the eccentricity of the binary companion, it is typically on the order of 
Gyr. Importantly, the Kozai mechanism takes place dominantly in single-planet systems, since mutual interaction 
between multiple planets usually suppresses the secular perturbation. Th e four single-planet s ystems in the regio n 
t > 9 Gyr in and e > 0.6 in Figurelll HD20782 (e = 0.92. Uones etaLll2006l). H D45350 (e = 0.78.' lMarcv et al.ll2005bh . 
HD222582 (e = 0.76. IVogt et al.l l2000f ) and HD3651 (e = 0.63. iFischer et al.l[2003h are also good candidates where 
secular eccentricity oscillations could have taken place. Since the Kozai mechanism can be effective in a fairly wide 
binary (a > 1000 AU), it is still possible that a low-mass companion for these systems has rem ained undetected. The 
results from ongoing searches for wide-orbit companions around extrasolar planetary systems (jMugrauer et all 120041 : 
iChauvirTe t al. 2006J) will help resolve these questions, particularly for old single-planet systems with very high orbital 
eccentricities. It should also be brought to attention that in Figure [ll?] there are at least a few single-planet systems in 
which large (~ 0.6) eccentricities have been excited within a few Gyr. Finding systems with very large eccentricities 
(e > 0.8) and young ages would be a challenge for the Kozai mechanism. 

6. SUMMARY 

We have calculated theoretical stellar parameters (M, r, R, M ce , R ce , log g, etc.) for 1074 stars from the SPOCS 
catalog. Using Bayesian analysis, we have adopted an appropriate choice of the a priori stellar parameter distributions 
and computed posterior PDFs for each parameter. We have provided several statistical summaries for each posterior 
PDF, such as the median, the mode, and various ranges of credible intervals, as well as flags indicating those cases for 
which a parameter is "poorly-determined" . The newly determined physical properties of the five sample stars proved 
to be consistent with the previous calculations by FRS99, but the high resolutions of our stellar evolution database 
provided stellar models with much smaller uncertainties (§ 15. ip . The complete list of the derived stellar parameters 
available in the electronic version of the paper are now ready to be used for various dynamical and formation studies of 
planetary systems. Also, the computed database of stellar evolution tracks will continue to be a useful tool for modeling 
extrasolar planetary systems. For example, the precise determinations of stellar radii applied to the large number of 
transit detections anticipated by OGLE or Kepler can provide important information for the interior structures of 
gas-giant planets. 

The uniformly analyzed physical properties of the large stellar sample provided several interesting relations between 
stellar parameters. The derived relations between convective zone mass and metallicity of known planetary systems 
seem to reject the atmospheric metal enrichment of planet-host stars caused by planet accretion. We found no 
significant evidence of stellar surface metallicity diluted in the large convective envelope of evolved stars, which was 
expected in the planet-pollution hypothesis. 

The sample size of known planetary systems in the SPOCS catalog still limits statistical analyses aiming to identify 
correlations between stellar and planetary parameters. The sample is mostly abundant with G-dwarfs due to their 
strong suitability for Doppler surveys. A greater number of K- and F-dwarfs and evolved subgiants are needed to 
increase the leverage for identifying correlations with stellar mass and radius. Nevertheless, some of the features 
observed in the scattered plots in Figure 1X711X91 might provide clues about the history of formation and dynamical 
evolution in planetary systems: 

(i) There is a lack of close-in planets (P < 100 days) detected around subgiant stars with masses greater than 1.5M ffi 
(Figur e 1X8 ]) . This may be an effect of small number statistics. However, it has been demonstrated by IRasio et al.l 
(1996) that the timescale for orbital decay due to tidal dissipation becomes progressively shorter as the planet-host 
star evolves into the subgiant phase. The five long-period planets around subgiants may indicate that many of the 
short-period planets suffer orbital decay and will eventually be engulfed in the atmosphere of the giant stars. Current 
planet searches around A - F stars and subgiant stars will help constrain the possibilities of such dynamical process. 

(ii) Many of the planets with extremely high eccentricities (e > 0.6) are discovered in old systems (Figure IX9|) . 
Most of them are in systems with only a single known planet, and two of them are considered very likely to have 
been affected by secular perturbations induced by a stellar binary companion. Secular orbital dynamics are likely to 
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be responsible for such systems, since any impulsive perturbation can disrupt the planetary orbits in much shorter 
timescales. It is of great interes t to see whether the co-proper motion sury eys in search for distant companions around 
known planetary systems fe.g.. iMugrauer et al.ll2004l ; IChauvin et al.ll2006f ) will confirm, or rule out the possibilities of 
stellar companions around these systems. 

This work was supported by NSF Grants AST-0206182 and AST-0507727 at Northwestern University. E.B.F. 
acknowledges the support of the Miller Institute for Basic Research. 
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